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ABSTRACT 


Three methods for weighting an exponential regression 
model to estimate discrete reliability growth were derived 
and tested. The first method systematically applies greater 
weight to test phases whose estimates have less variability. 
The second method similarly applies heavier weight to the 
most recent test phase estimate. The third method allows 
the user to choose the weighting scheme. 

These methods were evaluated against eight patterns of 
actual reliability by altering a previously developed Monte- 
Carlo simulation. Their performance was then compared to 
the unweighted exponential regression and Maximum Likelihood 
Estimate With Discounting (MLEWD) models. The second 
weighting method appears to perform the best under rather 
general constraints. Also, comparison is made of the least 
squares estimates for reliability growth using two different 
unbiased estimates for the negative of the natural logarithm 


of the failure rate. 
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I. INTRODUCTION 


Systems and equipment often go through several phases of 
development during which the reliability is assessed and the 
components or design are appropriately modified. From these 
early assessments, project managers need to make predictions 
regarding the ability of the system or equipment to meet 
reliability specifications by the contracted acceptance or 
delivery date. Some conventional estimators of reliability 
disregard data related to immature systems and use only that 
data obtained from the final product in verifying 
reliability. This practice is extremely inefficient, costly 
and occasionally, infeasible. Using all available test data 
would not only permit a manager to monitor the development 
process more successfully but at a substantial savings to 
all parties involved. 

Reliability growth models are useful in estimating the 
reliability of both immature and mature systems. Continuous 
and discrete versions of models exist. Models which employ 
attribute test data are discrete while those based on time 
to failure data are continuous. Only discrete models will 
be addressed in this thesis. Both of these versions make 
use of all available test data and are often used in 
conjunction with a technique known as failure discounting. 


This iS a process whereby an increasing fraction of a 


failure is repeatedly diminished as greater amounts of 
Subsequent data is accumulated without repeat of the same 
failure. 

Models evaluating data which have been adjusted using a 
failure discounting scheme have performed well in estimating 
actual reliability [Reft. b:pp. 852-55]. Relying on 
discounting techniques has its drawbacks however. In order 
to employ the straight percent discounting method, for 
example, the user must select two parameters, N and F. F is 
the fraction each failure will be discounted and N is the 
number of successive successes which must be achieved prior 
to applying the discounting method. No rules have been 
developed as to how these parameters should be selected and 
no guidance other than to "use good engineering judgment" 
{Ref. l:p. 50} has been offered on the subject. 

The objective of this paper is to provide the user of a 
reliability growth model with an alternative to failure 
discounting which will produce comparable results. To this 
end, the Monte-Carlo simulation developed by Captain James 
Drake and modified by Captain James Chandler was altered to 
include a weighted exponential regression model. This model 
can use one of three methods of weighting. 

Method one is based on the idea that the variances of 
the observables are not all equal and hence some are more 


reliable than others. To accommodate this fact, weights in 


this method are "“heavier" for observables with lower 
variances. Method two systematically gives more weight to 
the most recent test phase motivated by the idea that the 
most recent version of the system is likely to be the most 
reliable. Method three is an extension of this idea. In 
this method, the user is permitted to subjectively select 
how much weight he desires to give each phase. This might 
produce successful results if the user has an intimate 
knowledge of the development process. He could then 
possibly choose to give little weight, say, to a pnase he 
knew was conducted poorly or contained problems that he was 
certain were since corrected. Each of these methods was 
evaluated and compared to previously developed models using 
eight different actual reliability growth patterns. These 
comparisons are contained in Chapter VI. 

As a further modification to previous work which 
primarily explored the "test-fix-test" methodology, i.e. 
only one failure per phase of testing permitted, this thesis 
investigates the effect of using an alternative unbiased 
estimator of the exponential regression parameter of a 
"test-find-test" scenario. 

The following chapter will address previous work in the 
area of reliability growth models. The motivation for 
altering the exponential regression parameter is discussed 
in Chapter III and the weighted Regression Model and the 


derivation of the three methods of weighting are developed 


in’ Ghapter “Iv. A brief description of the models with 
which the weighted model will be compared is presented in 
Chapter V. The comparisons under varying patterns of actual 
growth is treated in Chapter VI. A summary, conclusions and 
recommendations for further study are discussed in Chapter 


ILO 


II. PREVIOUS WORK 

A. BACKGROUND 

This thesis is the third in a recent sequence on the 
subject of reliability growth. The initial work was 
completed by Captain James Drake [Ref. 1]. In his study, 
Captain Drake developed a Fortran program to evaluate the 
performance of three reliability growth models; namely, the 
maximum likelihood estimate model, the exponential 
regression model, and the weighted average model. These 
models were evaluated in conjunction with two failure 
discounting methods, the confidence limit (C.L.) method and 
the straight percent discounting method. The C.L. method is 
referred to as the Lloyd discounting method in Drake's 
thesis. In addition, Captain Drake varied the parameters of 
each of the discounting methods to evaluate their effect on 
the performance of these models. The "actual" or known 
reliability growth pattern against which these models were 
compared was generated by Monte-Carlo simulation. In this 
Simulation, the user inputs the reliability during the first 
phase of testing of each component which may possibly fail 
(i.e., Of each potential failure cause). The program then 
generates an "actual" growth pattern from these values for 
the remaining phases of testing. The user of this program 
relinquishes strict control of the actual reliability growth 


pattern beyond the first phase and, although manipulation of 


another parameter (specifically, "FRIMP") can produce the 
pattern in a general desired form, the simulation is 
incapable of producing a decreasing pattern. 

Captain James Chandler modified Captain Drake's program 
to permit the user complete control of the “acumem 
reliability growth pattern at each phase of development. In 
his version of the program, the user is required to input 
the reliability of each potential failure cause at each 
phase. Additionally, if more than one failure is permitted 
per phase of testing, the value must be replicated for each 
failure in that phase [Ref. 2]. To reduce the amount of 
required user input, this process was modified in this 
thesis so that the amount of information required per phase 
is the same regardless of the number of failures per phase 
permitted. The required information is contained in the 
User's Guide, Appendix A to this thesis. 

Captain Chandler also altered the C.L. failure 
discounting method to allow the user control of the discount 
interval. This was motivated by the results of the original 
thesis which indicated that employment of the C.L. method 
resulted in overly optimistic estimates regardless of the 
model used or the actual reliability pattern generated. 
This was attributed to the fact that the original C.L. 
method diminishes the weight of a failure after each 


following success. 


B. RESULTS OF PREVIOUS WORK 

The weighted average reliability growth model was 
eliminated as a potential model in the early stages of the 
original analysis. This was due to its consistent 
overestimation of actual reliability which was only 
amplified by employment of either of the failure discounting 
methods. This model was discarded in subsequent work as 
well. 

The remaining two models accurately tracked a wide range 
of reliability growth patterns. The Maximum Likelihood 
Estimate with Discounting model, (MLEWD), generally “tended 
to underestimate actual reliability in early phases and 
slowly converge to the actual value with increased test 
Gata" [Ref. l:p. 48]. This model also exhibited smallest 
variance of all models evaluated. The negative aspect of 
the MLEWD model is that the choice of discounting parameters 
is critical to successful implementation and guidance in 
their selection is non-existent. The potentially drastic 
effects of different parameter choices is discussed in 
Chapter V. 

The exponential regression model also performed well 
against various actual reliability growth patterns and was 
far less sensitive to the choice of discounting parameters 
than the MLEWD model. The difficulty with this model is 
that it is highly variable in the early phases of testing 


although it generally stabilized after four phases. This 


was considered acceptable since most development processes 
allow for more extensive testing. 

The C.L. discounting method was not recommended for 
employment with either model as its use only occasionally 
produced results comparable to the straight percent 
discounting method without possessing its corresponding 
flexibility. 

Captain Chandler evaluated these two models along with 
the standard single phase maximum likelihood estimate 
against eight reliability patterns. His results indicated 
that both the MLEWD and the exponential regression model 
were Superior in all respects to the standard estimate. 
Additionally he found that the modification made to the C.L. 
discounting method resulted in instances of Superior 
performance to the original C.L. method (although not to the 
straight percent discounting method). The remainder of his 
work substantiated the conclusions drawn by Captain Drake 
and included additional observations regarding the case of 
declining reliability. 

Both authors conducted all simulations for the test-fix- 


test scenario. 


IIIT. MOTIVATION FOR IMPLEMENTING ¥;. 


The exponential regression reliability growth model was 
developed by H. Chernoff and W. M. Woods. Its derivation 
is fully detailed in Ref. 3 and Chapter IV. In this 
model the reliability after the k*® change is modeled by 
R,.=1-e" (@+6K) | For the purposes of this chapter it is 
sufficient to understand that the model estimates the 


reliability in Phase k as: 


— ] — on Ort Bid) 


(3.1) 


where Ry, is the reliability estimate, and the estimates @, 
and By for a and 8 at the conclusion of testing in the qth 
phase are obtained using linear regression methods and an 
unbiased estimator for or + By) - 


The unbiased estimator is: 


— fle df24+ 1/3 +...0/(X%, —1) ior Xj, =e 
Yk = 0 : Mine Kike 


(3:2) 
where Xa is the number of trials between the (3-1) St 
failure and the 3th failure (including the jth failure) in 


the «th phase. 





* The symbol B and B denote th 
@€ same constant th 
this thesis. B is used for 8 within equations. - 


Let F, = the number of falluressimgere qth phase. Since 


Xq ke 2 ak eae , are independent random variables then 


Y, S@4peaioe YER (3.3) 


is also unbiased. 

Previous work incorporates this estimator of (a + 8,k) 
however an unbiased estimator which has minimum variance of 
all estimators has since been developed by W. M. Woods; 


this estimator is as follows: 


Y,.* = Fig t U(Fy +1) + U(Xp, — 1) (3.4) 


Replacing Yy with Yj} in the exponential regression model 
provides a more accurate estimate of reliability. 

Table 1 is the result of performing ten thousand 
replications of one phase negative binomial (Pascal) test 
data with actual reliability R and comparing the model to 
this value of R employing Yy and again using Y),- This table 
demonstrates that while reliability predictions using 
Yj and Ge are both conservative, those using Y} more closely 
estimate actual reliability, have smaller variances and a 
smaller mean Square error. 


The parameters used in this simulation are as follows: 


F, = The number of failures in phase k was fixed atuomaas 
all cases 
A ee number of trials to the ie failure. 
k 


10 


The number of replications = 10,000 


Z 
ll 


nO 
ll 


The number of phases = 1 


Yy and Y3 are as in equations 3.3 and 3.4 respectively. 


TABLE 1. COMPARISON OF Y,; VERSUS Y, AS THE ESTIMATOR OF A, 


Actual 
Reliability 


0.9269 
starsat [tease | _3c6s3_[ ozs] 02367 |_ 02a 





Though these results prove that Yy is superior to Yy as 
an estimator of the coefficient of the exponential term, the 
impact of this improvement is disappointingly trivial. 
Figures la - lh depict the results of employing y* and Y to 
the exponential regression model for actual reliability 
patterns or As one can see from both the Table and the 
figures, the effect of this change takes place in the third 
decimal value. In practice, changes of this magnitude are 
hardly useful. Still, since — Y;. is an improvement over 
Y, it is retained as the method of estimating (a and 8k) 
throughout the remainder of this thesis. 


Note: In all figures, A=actual; *=y* and - = ¥ 
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IV. THE WEIGHTED EXPONENTIAL REGRESSION MODEL 


A. THE EXPONENTIAL REGRESSION MODEL 

The exponential regression model obtains sequentially 
updated estimates R, of the reliability R, which denotes 
true reliability after the qth phase. The basic model for 
R, is: 


R,=1—e7 4” (4.1) 


In the exponential regression model, linear regression is 


used to estimate A,= (@+8k) by A, where 
So on + Byk (4.2) 


The exponential regression estimate of reliability is then: 


Ry = 1 —e7 Het BM k= 1,2 (4.3) 


ge ae 


suits Osseo: Byk is estimated sequentially at each phase. This 
results in a model which is capable of tracking changing 
reliability. In order to estimate the parameters a, and By 
one must first calculate the unbiased estimator, Y}- As 
discussed in Chapter III, Yy 1s the unbiased estimator with 
minimum variance and has been implemented in this version of 
the model. 
er: 


Fr »= the number of failures any the qth phase 


ras 
ry 
ll 


The total number, of trials to and including @ chemaes 
k failure in the kth phase. 


iG 


Then: 


Y,* = aa WF +1) + A/(Fy + 2) + ..1/(Xp, -1) epee 
jor aE (4.4) 


As an example, if the testing in the fifth phase 
continues until three failures occur and the third failure 


occurs on the eighth trial then: 
Ye = 1/3 + 1/4 + 1/5 + 1/6 + 1/7 = 1.092857 


The least squares estimates for 8, anda, are : 


k 
yc —k)x ¥*,] 
= 





Bo = 
K k : (4.5) 
a: 
Fa 
Co, = ‘~— B,k (4.6) 


where k=(1 +24+34..4)k 


¥,=(¥*, + Y*, + Y*,4..¥*)/k 


Replacing the unknown parameters, a, and 8, with their 


estimates @, and By in equation (4.3) yields: 


A as — (3 + B,k) 
R, = or for a 
e for k=] (4.7) 
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B. THE WEIGHTED EXPONENTIAL REGRESSION MODEL 

The weighted exponential regression model is identical 
to the unweighted exponential regression model above with 
the following exceptions: 


Let: 


w; = the weight applied to the estimate in@phasouam 
Replace all occurrences of k in the calculation of the 
unweighted model with: 


k 
>= 9 


i=] 


k, =———_ 


(4.8) 
S 


=| 


and replace ve everywhere with: 


k 
’ "Xe 
a ; 
a —_ i=] 


kw — n 
>, (4.9) 


Making these substitutions, the estimates of a and 8 are 


now ; 
k 
> [G~&) x ¥f x wi 
By = (4.10) 
» li —k,) x w, 
=! 
orw = Yew 7 Buwkw (4.11) 


A _ fi e7 Gan + Bind | for k> I 
Raw = [! Ser S| Se for k=1as before (4.12) 
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C. EXAMPLE OF THE WEIGHTED EXPONENTIAL REGRESSION MODEL 

The following example is offered as a means of 
Clarifying the application of this model. 

Suppose that a system undergoes ten phases of testing 
and a phase is terminated upon occurrence of the third 
failure in each phase (i.e., F,=3, Fo=3..-- Fyi9=3). Suppose 
further that the weights given to each phase are as follows 


(See Table 1A): 


TABLE 1A. EXAMPLE OF WEIGHT DATA 










| Phase | 1 | 2 | 3 | 4 | 5 | 6 | 7 {8 | 9 | 10 | 


These weights have been arbitrarily selected for the 
purposes of this example. Actual calculation and/or 
selection of weights will be discussed in section D of this 
chapter. 

Tables 2 and 2A represent the test data collected on 


this fictitious system and the resulting calculations based 


on this data: 


i, 


TABLE 2. EXAMPLE PARAMETER CALCULATIONS BY PHASE 
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and reliablilitv estimates, R, = 1 — e7 Gu ~ Bra! 


TABLE 2A. EXAMPLE RELIABI EGY Ese 










Phase (kj 
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= not computed in the first phase. See equation (4.7). 
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D. METHODS OF WEIGHTING 


1. Method One. 


In the unweighted exponential regression model the 
estimate of 8 is the Gauss-Markov estimate when the variance 
in each phase is equal. When the variance is not equal from 
phase to phase "the variance of 8 is unnecessarily large." 
(Ref. 5]. This variance in parameter 8 accounts for the 
demonstrated higher reliability estimate variance evidenced 
in previous work [Ref. 1]. In order to diminish this 
variability each phase is systematically accorded a weight 


determined by the estimate of the variance of the phase: 


Ww, = the weight allocated in phase k = Sa eee a (4.13) 


using weights derived in this manner, %,,, in equation (4.9) 


becomes: 


08 (Geshe x ta 


Baw = (4.14) 


, 
> Le — Ky)? x 1/673 
= 


which is unbiased and has minimum variance among all 
unbiased estimators of 8 [Ref. 5]. 
Figure 2 demonstrates the effect of weighting the 


exponential regression model by this method. 


Ze 


As one can see, the weighted model more closely 
estimates actual reliability, particularly after the third 
Phase of testing. Chapter VI contains more detailed 


comparisons of models. 


A=ACTUAL, R=REGRESSION, 1=METH.1 


1.0 


0.8 


RELIABILITY 
0.6 


0.4 





2 4 6 8 10 
PHASE (1 FAILURE PER PHASE) 


Figure 2. The Weighted Exponential Regression Model 
in Conjunction with Method One Weights. 
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2. Method Two 

The second method of weighting applied to the 
exponential regression model was less theoretically derived. 
This method systematically allots more weight to the most 
recent phase. The rationale for the implementation of such 
a method is that since the aim of the development process is 
to improve the system, the reliability in the latter phases 
will be greater than in earlier phases. While this is not 
necessarily the case, the method none the less has intuitive 
appeal. 

The weights used in this second method are 


calculated for each phase as follows: 





(4.15) 


To see why this expression creates weights which 
increase with phase, recall that an estimate of the variance 


of the negative binomial distributed data is: 


A F,) x (p, 
on = an estimate of variance in phasek = ae (4.16) 
(Pre 
where: 
F. = # of failures in phase k, as before, and 
Py = The maximum likelihood estimate of the 


reliability in phase k. 
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Assume, aS in our example in section C, that the 
number of failures in each phase is constant, then the term 
F. may be ignored. One can readily see that as Py 
increases, ae increases and therefore, w, increases. 

Figure 3 demonstrates the effect of weighting the 
exponential regression model by this method. The model 
weighted in this manner more accurately tracks actual 
reliability than both the unweighted and method one weighted 


model. More extensive comparisons are drawn in chapter VI. 


A=ACTUAL, R=REGRESSION, 1=METH.1,2=METH.2 


RELIABILITY 





4 6 8 
PHASE (1 FAILURE PER PHASE) 


Figure 3. The Weighted Exponential Regression Model 
in Conjunction with Method Two Weights 
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3. Method Three 

The third and final method of weighting the 
exponential regression model has no roots in theoretical 
mathematics. This "method" is simply to let the user select 
the amount of weight he desires to assign to each phase. 
His only constraint is that the sum of these weights over 
all phases equals unity. Selecting the weights for this 
method is currently analogous to specifying the parameters 
for use in failure discounting. Unless one possesses 
extensive knowledge and intimate familiarity with the system 
being evaluated as well as its developmental history, use of 
this method is not advisable. 

To demonstrate the radical behavior of the model 
when this method of weighting is employed, five different 
cases of weights were evaluated. The weights for each case 
are listed in Table 3. 

Figure 4 is an example of the behavior of the model 
when the weights listed as case 3 are employed. As one can 
see, this selection of weights resulted in an extremely poor 
performance of the model against the pattern of actual 
Beliability. As with the previous two methods, further 


comparisons are contained in Chapter VI. 
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TABLE 3. SAMPLE METHOD THREE WE2ZGHTs 







1 | oraassria | oiazssmia | 03 | 333333933. | 0005 
4 ovazasrva | ova [05 sss [oes 








A=ACTUAL, R=REGRESSION, 1=METH.1,2=METH.2,3=METH.3 
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Figure 4. The Weighted Exponential Regression Model 
in Conjunction with Method Three Weights. 
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As a comparison, the weights for each of the methods 


depicted in Figure 4 are listed in Table 4. 


TABLE 4. COMPARISON OF METHODS ONE, TWO AND THREE WEIGHTS 
OVER A TEN PHASE SIMULATION. 










ee Ea a a a a a 
Method | , P 
O68? P89 HOSTS OOI8? 
Nlethod 2 $55 
Weights 1.0 | S00; o> S30E74 1297099 1279019 2192977 | .4OS804d3 | 2705978 | 695464 
Nlethod 3 
4 c S § 2 










The values listed are the amount of weight given to 
that particular phase when it is the current phase. Thus, 
for example, in Phase 9, method 1 allots .0031794 to the 9th 
phase and 1- .0031794 = .9968206 to the previous eight 
phases. Method 2 allots a much greater weight to phase 9, 
Mee., -705975, and only 1-.705975 = .294025 to the previous 


eight phases. 


Pad 


V. THE MAXIMUM LIKELIHOOD ESTIMATE WITH FAILURE 
DISCOUNTING IN THE TEST-FIND-TEST SCENARIO 


A. DESCRIPTION 
The conventional maximum likelihood estimate of 


reliability is: 


R= # of successes/# of trials (Sa 
In order to use the estimator however, "...constant 
reliability, R, is required for each trial. Because the 


reliability at each phase R, may not be constant with k; 
only the test data from the phase of interest may be used to 
estimate reliability" [Ref. l:p. 22]. An excellent detailed 
description of failure discounting as applied to the test- 
fix-test scenario is contained in Ref.1l:pp.13-21. This 
section primarily addresses the process when employed in 
conjunction with the test-find-test scenario. 

The straight percent discounting method attempts to 
reflect improved system reliability by removing a fraction, 
F, of a failure’s weight at an interval of every N trials. 
This sequence of N successful trials must occur in a follow- 
on phase of testing, after the components which were "found" 
to be causes of failure are modified or repaired. The idea 
here is that accumulation of subsequent testing without 
repeat of a failure for the same cause increases the 
confidence that the failure cause has been removed; 


therefore its weight in further estimations should be 
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diminished. When using this method, the data is first 
adjusted and the model is applied to the modified data. The 
Success of the model is wholly dependent on the correct 
adjustment of the data and hence the parameters chosen in 
calculating these adjustments. The equation used to compute 
the current value of an adjusted failure is: 

Adjusted failure (;) = (1-F) int (M/N) (5.2) 
where M is the number of successful sequential trials in 
follow-on phases for failure cause j} and F and N are as 
above. 

Pee EXAMPLE 

As an illustration of the method, consider the results 
of the following example. Suppose, as in our previous 
example, testing in each phase continues until the 
occurrence of the third failure. Further suppose the 


results of two phases of testing are as in Table 5. 


TABLE 5. MLEWD EXAMPLE DATA 


# Failures # Trials to 3rd Failure 





PS, 


Then applying equation (5.1), the maximum Pekelihees 


estimates without failure discounting are: 


Rye) oe R> =) C/ eee eee (5035) 


As one can see, the information obtained from the first 
five trials is disregarded in the calculation of the 
estimate of the reliability in the second phase. Continuing 
in this manner, the data in the first eleven trials will be 
ignored in computation of the third phase reliability. This 
practice is extremely inefficient by any standards. 

In the test-fix-test scenario, [Ref.l:pp. 15-16] it is 
assumed that a design "fix" is implemented after each 
failure; ina test-find-test scenario, no modifications are 
assumed until a given number of failures is observed. The 
causes of the failures are merely noted for subsequent 
action: Since there is no justification for applying the 
discounting method until follow-on phases prove a cause 
corrected, different values of weights are assigned to the 
Same sequence of successes and failures in these two 
scenarios. Tables 6 and 7 illustrate the application of the 
discounting method to both situations. In the calculations, 
the parameters N=3 and F=.25 are arbitrarily chosen. 
Unfortunately, in practice, they are often likewise 


selected. 
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TABLE 6 


TEST-FIX-TEST DATA WITH DISCOUNTING 


FAILURE CAUSE 


Phase Trial A B C D 
al 1 S O S O S 0 F 1 
Z 2 S O F 1 S 0 S 1 
3 3 S 0 S 1 S O S 1 
3 4 S 0 S 1 S O S Re gs) 
3 5 F 1 S </o 5 O S 75 
4 6 S 1 S aro 0 S hs 
4 7 S 1 S wot ce: 1 S no 5 
5 8 S Be) eS 5625)| S it F 1 
6 9 S ; Slee 5625] S 1 S 1 
6 10 S B/5| F JL S 2/39 2s ik 
i dag S 5625/| S 1 S e7 ON S B73 
Netation: For each failure cause the attribute 
S = Success F = Failure 


is listed along with the corresponding weight of that 
attribute. 
TABLE 7 
TEST-FIND-TEST DATA WITH DISCOUNTING 
FAILURE CAUSE 


Phase Trial A B Cc D 


ono 
Sa Oo ee 2 © 


NNNNNNAFRPRE FB 
FOUDNAUAWNE 
ANNNNNWHANnNHnHN 
NAMNNNNNHNNHnHWNHN 
PREYINEP RIP RPRPRO 
NANNNMANHNHNHAN 
ANN WWHNnNANNNnANnNN 
PRR RP RP RIP RP RPE e 


ioe 


Cul 


Testing in the Test-Find-Test scenario terminateswatem 
the 3rd failure. The straight percent ~discouneum, 
parameters are: 

N = Discount Interval = 3 


1g 


Fraction the failure is reduced = .25 


Using the test-fix-test discounted data and MLEWD model 
produces the results shown in Table 8 upon completion of the 
eleventh trial. 


TABLE 8. RESULTS OF THE MLEWD MODEL APPLIED TO THE 
TEST=FIA=-TEsT pA 


PHASE = 
ULC lua A CAUSE ADJ. FAILURE ADJ. TRIALS 





Dt IES. ls 0005. 


Adj. trials = (the observed # of trials since the previous failure); 
( adjusted failure value). 


and 


. © ADIUSTED TRIALS) =(@OTAL HOP Buus 
Rg = eS = 909 (5.4) 
() ADJUSTED TRIALS) 


3iZ 


As one can see, this estimate is higher than that 
predicted after the eleventh trial (second phase; eq. 523) 
using the MLE without discounting. Since it takes into 
account 5 additional data Holmes gee iS considered to be a 
superior estimate...provided the parameters were correctly 
chosen. 

The results of the test-find-test-scenario after the 
eleventh trial are as shown in Table 9. 


TABLE 9. RESULTS OF THE MLEWD MODEL APPLIED TO THE 
TEST-FIND-TEST DATA 


PHASE CAUSE ADJ. FAILURE ADJ. TRIALS 
po 1.0 = 1.0000 
po Bt. = 1.0000 
pA 5625 35625 = 5.3333 


_ eee oe 
C 21.0 = 2.0000 
3.1.0 = 3.0000 





mes. TRIALS = 13.3333. 


A 


R, = .5499 


As a comparison, the reliability after the second phase 
of testing in the test-fix-test scenario would be zero since 
there were two failures in two trials. 

The results seem to indicate that the estimate of 


reliability of a system or equipment will grow more Gunek ly 


S13) 


appeal as it seems likely that if failure causes are 
corrected immediately they are else like to be responsible 
for further failure. These results additionally demonstrate 
that the use of discounting has the desirable characteristic 
of producing higher estimates of reliability than the 
conservatively biased MLE without discounting [Ref. 3:p. 


34). 
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VI. MODEL COMPARISONS 


A. INTRODUCTION 

In sections B through F of this chapter, the performance 
of the unweighted exponential regression model will be 
compared to that of the weighted model used in conjunction 
with weighting methods one and two. The employment of 
method three will be discussed in Section G. The 
performance of the MLEWD will additionally be discussed 
where appropriate. For the purposes of this evaluation 
patterns VI, VII and VIII will be discussed aS a group in 
the category, Constant Reliability. Patterns IV and V will 
be Similarly grouped under the title Rapid Reliability 
Growth, while pattern I, Convex Reliability, pattern II, 
Decreasing Reliability and pattern III, Intermittent 
Reliability will be evaluated individually. All discussion 
is in regard to test-fix-test scenarios to facilitate 


reference to previous work. 


B. CONSTANT RELIABILITY PATTERNS 

Figure 10 [Ref. 2:p. 32] illustrates the performance of 
the MLEWD model against representative constant reliability 
pattern VI. No failure discounting was implemented in this 
case. The MLEWD model performed better than either the 
weighted or unweighted regression models in all three 


constant reliability patterns, although its superiority 


35 


diminished with lower levels of constant reliability. This 
was expected since, as indicated in Chapter V, the 
underlying assumption of the MLE is one of constant 
reliability. 

Figure 5 is an illustration of the unweighted and 
weighted exponential regression models also against constant 
reliability pattern VI. Regardless of the level of 
constancy, the unweighted model consistently underestimated 
reliability while the weighted model overestimated it. in 
all patterns, the model weighted using method 2 (i.e., more 
weight to the most recent phase) produced better results 
than when method 1 (i.e., more weight to less variable 


phases) was employed. 


C. RAPID RELIABILITY GROWTH 

Figure 6 depicts the performance of the models when 
applied to Pattern V. This pattern is one in which "actual" 
reliability achieves its greatest value (approximately .9) 
early in the testing process (phase 3) and remains there 
throughout the remainder of the testing. Pattern IV is, in 
all ways, Similar to pattern V with the difference being 
that it achieves the value of approximately .99 in phase 3. 
Both versions of the weighted model and the unweighted model 


underestimated reliability during its growth phases although 
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Figure 5. Weighting Methods One and Two and the Constant 
Reliability Pattern 
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Figure 6. Weighting Methods One and Two and the Rapid 
Reliability Growth Pattern 
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the weighted models more accurately tracked the process. As 
with constant reliability, the model weighted by Method 2 
outperformed the model when weighting Method 1 was employed. 
Also similar to constant reliability performance results, 
both versions of the weighted model produced higher 
estimates than the unweighted model over all phases. Given 
this underlying pattern, if one is concerned with tracking 
the initial phases of testing, the weighted model in 
conjunction with Method 2 is recommended as it best captures 
the growth portion of the development process. The 
unweighted model in this case seems to best capture the 
constant reliability portion of the development although its 
superiority is minimal. Figures 13 and 14 [Ref. 2:pp. 36- 
37] illustrate the performance of the MLE both with and 
without discounting. In both cases, the exponential 


regression model outperforms the MLE over all phases. 


De CONVEX RELIABILITY GROWTH 

This pattern illustrated the most dramatic support for 
application of weighting to the exponential regression 
model. Figure 7 depicts the substantial improvement in 
reliability estimation when weighted Method 2 is employed. 
Use of Method 1 also produced superior results to the 
unweighted model although these were not as dramatic, 
particularly in the latter phases of testing. Both versions 


of the weighted model seemed to capture the pattern of 
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CONVEX RELIABILITY GROWTH PATTERN 
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Figure 7. Weighting Methods One and Two and the Convex 


Reliability Growth Pattern 


growth as early as phase 2. No earlier improvements can be 
anticipated as all three versions of the regression model 


require at least two phase of testing to perform a linear 


regression. 


Figure 23 [Ref. 2:p. 46] demonstrates the poor 


performance of the MLE when no discounting is applied. 


SN 


These resultS are not surprising as the underlying 
assumption of constant reliability in the case of the MLE is 
violated here. A dramatic improvement as seen in Figure 24 
[Ref. 2:p. 47] when discount parameters, F=.75 and N=6 are 
applied to the data. The parameters were selected after a 
good deal of trial and error simulation, however, and no 
explanation as to why this choice of parameters was 
effective can currently be offered. Since the luxury of 
simulation is not often available in practice, and the 
weighted model in conjunction with Method 2 outperforms the 
other models, it is recommended for use while the underlying 
reliability pattern is suspected to be convexly increasing. 
E. DECREASING RELIABILITY 

Figure 8 indicates that the same general previous 
observations can be made with regard to the performance of 
the models in the case of decreasing reliability. The 
unweighted regression model produced estimates which are 
everywhere below those of the weighted model regardless of 
the weighting method employed. Additionally, the weighted 
models seem to more accurately capture the changing growth 
pattern. The unweighted model is most accurate when actual 
reliability decreases to a low value however this is 
primarily due to the fact that the unweighted model 


consistently tends to underestimate actual reliability and 
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DECREASING RELIABILITY GROWTH PATTERN 
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Figure 8. Weighting Methods One and Two and the 
Decreasing Reliability Growth Pattern 
is less responsive to changes in the actual growth pattern. 
Therefore when actual reliability dipped and the already low 
unweighted model did not eS euel. the result was a fairly 
accurate estimate. As in previous patterns, the model 


weighted by giving most weight to the most recent phase 
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estimate produces the most responsive results. For this 
reason, this model would be recommended in this case. 
Figures 17 and 18 [Ref. 2:pp. 40-41] demonstrate the 
inability of the MLE model to capture a decreasing 
reliability pattern whether failure discounting is employed 


Or NOt. 


Ee INTERMITTENT RELIABILITY GROWTH 

The Intermittent Reliability Growth pattern (Figure 9) 
can be described as one in which reliability grows fairly 
rapidly, remains constant for several phases and then 
resumes growth. Such a pattern would occur if intended 
improvements to a developing system were ineffective for a 
period of time before discovery of beneficial change was 
implemented. As in all previous patterns, the unweighted 
model produced more conservative estimates than both 
versions of the weighted model. Also as before, the model 
weighted by Method 2 was most responsive to changes in the 
reliability pattern. This is particularly evident in phases 
59 and 6 where it dips below the Method 1 estimate in an 
attempt to recognize the period of constant reliability. 
All three versions of the exponential regression model 
outperformed the MLE model regardless of whether failure 
discounting was employed. Figures 20 and 21 [Ref.2 :pp. 43- 


44] reflect this. 
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INTERMITTENT RELIABILITY GROWTH PATTERN 
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Figure 9. Weighting Methods One and Two and the 
Intermittent Reliability Growth Pattern 


Ge THE WEIGHTED MODEL IN CONJUNCTION WITH METHOD 3 
WEIGHTING 


Weighting Method 3 is the application of user selected 
weights to the exponential regression model. Figures 10 - 
14 corresponding to representative examples of constant, 
rapidly increasing, convex, decreasing, and intermittent 


reliability respectively, demonstrate the sensitivity of the 
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model to the choice of weights. In every case the weighted 
regression model identically tracks the unweighted model 
until the user influences the process by applying a weight 
to a phase. For example, when the case 1 weights are 
employed, the user has given equal weight to the first seven 
phases of testing. This is equivalent to not weighting the 
model at all. Therefore, one would expect the unweighted 
model to produce identical results to this version of the 
weighted model. That this is in fact occurring, is most 
Clearly seen in Figure 10. Tracking the curve labeled "1", 
(corresponding to case 1, Weights Table 3), one can see it 
is identical to the unweighted growth curve until Phase 7 
after which it abruptly increases. Similar departures from 
the unweighted curve are seen for the other four cases of 
user selected weights. These results are welcome in that 
they substantiate the programming of the weighting process. 
They are disappointing however, in that, as with the MLEWD 
model, the response of the exponential regression model 
seems to be very sensitive to the choice of weights. 
Unlike the MLEWD model though, there are two systematic 
means of selecting weights which consistently produce 
accurate results (i.e., methods 1 and 2). It is recommended 
that these two weighting methods be employed unless the user 
is intimately familiar with the development process and has 
the unique ability to reflect his knowledge in the choice of 


weights. 
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Figure 10. Method Three Weights and the Constant 
Relivo. Groowe lerat fern 
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Figure 11. Method Three Weights and the Rapid 
Reliability Growth Pattern 
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CONVEX RELIABILITY GROWTH PATTERN 
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Figure 12. Method Three Weights and the Convex 
Reliability Growth Pattern 
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Figure 13. Method Three Weights and the Decreasing 
Reliability Growth Pattern 
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VII. SUMMARY, CONCLUSIONS, RECOMMENDATIONS 


A. SUMMARY 

Captain Drake created a FORTRAN program which simulated 
actual reliability patterns and evaluated the performance of 
three reliability growth models against these patterns. His 
work indicated that the MLEWD model was effective against 
all patterns when used in conjunction with the correct 
discounting parameters. He also found that the exponential 
regression model likewise tracked all actual patterns and 
waS relatively insensitive to the discounting parameters. 
This model, was however, more variable than the MLEWD model. 

Captain Chandler modified Captain Drake's program to 
produce several additional actual reliability patterns and 
Subsequently evaluated both models against these patterns. 
His results demonstrated that the MLEWD model was incapable 
of tracking declining reliability and was only Superior to 
the exponential regression model when actual reliability was 
constant; a phenomenon which, by design, rarely occurs in 
the development process. 

The objective of this thesis was to improve the 
exponential regression model by applying weights to the 
model at each phase. The intent was to produce an accurate 
means of estimating changing reliability without burdening 


the user with selecting discounting parameters. 
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Three methods of weighting were derived. ine -first 
method gives low weight to the most variable estimates of 
reliability, while the second weights most recent results 
more heavily. With both these methods, the FORTRAN program 
systematically computes these weights at each phase. [In 
method 3, the user is required to input values of weights 
for each phase of testing. Use of this method produced 
results which were very sensitive to the weights chosen and 
presented the user with the dilemma of selecting critical 
input parameters without any guidance. Since this is the 
very situation which the employment of weights was derived 


to avoid, this method is not recommended. 


B. CONCLUSIONS 

Regardless of the underlying actual reliability pattern 
or the method of weighting, the weighted exponential 
regression model produced higher estimates of reliability 
than the unweighted model. Since the unweighted model 
generally underestimates actual reliability this is 
considered to be a desirable characteristic. Both the 
weighted and unweighted model produced results comparable or 
Superior to the MLEWD model against all variations of the 
actual reliability pattern with the exception of the case of 
constant reliability. In this case, the MLEWD would be 


recommended. 
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In all other cases, the regression model weighted by 
Method 2 seems to be most responsive to changes in actual 
reliability. Additionally, it produced more accurate 
estimates than the model weighted by either Methods 1 or 3. 
In the case of convex reliability Method 2 weights are a 
clear choice over the unweighted model as well. When a 
period of decreasing reliability is suspected in the 
development process, the Method 2 weights are also 
recommended unless one is concerned with the magnitude of 
the decrease as opposed to the trend of the pattern. In 
this case the more conservative unweighted model is 
recommended. When estimating the reliability of the 
intermittent and rapidly increasing patterns, the Method 2 
weighted model most accurately captured the growth phases of 
reliability but were less accurate than the unweighted model 
in the final phase of testing although the differences in 
estimates here were minimal. Based on these conclusions, it 
is recommended that if a user has little or no Knowledge of 
the actual reliability growth of a developing system or 
equipment, he would be wise to select the exponential 
regression model in conjunction with Method 2 weights to 
estimate the reliability growth curve. 

As a final conclusion to this study, implementing the 
unbiased estimator with minimum variance, Yi in all versions 
of the exponential regression model resulted in an improved 


reliability estimate over the models in conjunction with Yue 
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RECOMMENDATIONS FOR FURTHER STUDY 


The following are recommendations for further study. 


Since a large segment of the commercial world collects 
its reliability data in terms of mean time between 
failures (MTBF's) it is highly recommended that a 
simulation to analyze continuous data be derived. 


Presently, straight percent failure discounting has 
shown some promising results when the discounting 
parameters are correctly chosen. Development of 
guidelines for selecting these parameters would 
certainly be of value in the employment of the MLEWD 
model. A difficulty in studying these parameters is 
that they cannot be altered from phase to phase or 
varied with failure cause in the current version of 
the program. 


As with failure discounting, the employment of user 
selected weights might prove beneficial if guidelines 
could be developed to aid in their selection. A more 
extensive analysis of this weighting method is 
recommended. 


ok 


APPENDIX A 


USER'S GUIDE TO DISCRETE RELIABILITY GROWTH (DRG) 


jae Introduction 
2. The DRG Exec File 


3. The Input Data File 
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APPENDIX A: USER'S GUIDE TO DISCRETE 
RELIABILITY GROWTH (DRG) 
fae LAtroduction 

In order to use the Fortran program, DRG, the user must 
possess three files: 

Te aenpue daca AL 
2. DRG Fortran Al 
3. DRG Exec Al 

A sample of each of these files is contained in Appendix 
B along with sample output. The input file and the exec 
file may be tailored to the user's needs. In its current 
form, the exec file produces a large degree of intermediate 
calculations for both the DRG Fortran program and Captain 
Chandler's program, JIMC Fortran Al. A copy of this program 
is also contained in Appendix B. These calculations may not 
be of interest to the user and may be eliminated with no 
detrimental effect to the program. A detailed explanation 
of each file definition is contained in section 2 of this 
appendix. 

The difference between the two versions of the Program 
is that JIMC Fortrans retains the ability to employ 
the weighted average estimate model and uses Y vice y" in 
the calculations of all versions of the exponential 
regression model. (see chapter III for a discussion of the 


motivation for altering this parameter which is discussed in 
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Section 3 para eZ. In all other ways, the input files for 
these two programs are identical. 
2. The DRG Exec File 

This file contains all the requisite file definitions 
and commands to run either DRG Fortran Al or JIMC Fortran 
Al. To indicate which of these programs is desired, the 
user must alter the third line of the file so that it reads, 
&FN=DRG or &FN=JIM respectively. Once this has been done, 
the user simply types DRG or JIMC while in CMS to execute 
the program. 

Of the remaining lines, only those beginning with the 
word FILEDEF should be altered or eliminated. All other 
lines contain commands pertaining to the execution of the 
program. The first FILEDEF line currently reads FILEDEF 10 
DISK INPUT1 DATA Al. The "1" in the filename indicates that 
this data file produces actual reliability growth pattern 1 
described in this thesis. Eight patterns and input files 
have been established for each program. These input files 
may be duplicated as described in Section 3 using the 
corresponding pattern parameters in Ref. 2 if the user does 
not have them readily available. This FILEDEF may not be 
eliminated from the exec file, however the "1" may be 
changed to any integer up to and including "8" for use with 
DRG Fortran Al. The integers 9 through 16 correspond to the 
Same patterns when executing JIMC Fortran Al. This 


numbering convention is adhered to with the remaining 
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FILEDEFs. Since a filename is only permitted to be six 
characters in length, whenever necessary, it is truncated to 
include the integers discussed. For example, input 1% is an 
unacceptable filename therefore the "T" is dropped from the 
word "input" to allow "inpul@." 

The next two FILEDEFs have the filenames Al and AQ 
respectively. The first corresponds to the “actual” 
reliability pattern (pattern 1) produced by the DRG program 
while the second corresponds to the same pattern when JIMC 
Fortran Al is run. These files should contain identical 
eucput. Both are kept for convenience so that if for 
instance, the user runs DRG with pattern 1 input and then 
JIMC with pattern 2 input he will still retain the initial 
"actual" reliability file, i.e., it will not be overwritten. 
These two files were both devised solely for plotting 
purposes and may be eliminated with no disruption to the 
program. The information contained in these files is 
duplicated in PRELIAB and JRELIAB Listing respectively. 

PRELIAB LISTING Al contains the primary comprehensive 
output for the DRG Fortran program while JRELIAB contains 
Similar information for JIMC Fortran Al. These files 
consist of a model parameter summary, a comparison of the 
mean predicted reliability for each model to the actual, the 
estimate standard deviation and a 95% confidence interval 


for each model. Also contained in these files is a 


ONS, 


recapitulation of results so that each model may be compared 
to all others. These FILEDEFs should not be eliminated from 
the exec file. 

JThesis and PThesis are the filenames corresponding to 
the intermediate summary output produced by JIMC and DRG 
Fortran respectively. These files contain a phase by phase 
listing of failures and failure causes for up to 5 
simulations. This output may be "turned on (or off)" by 
placing a 1 (or @)" in the appropriate line of the input 
file. See section 3 para 12 for details. These FILEDEFs 
should not be eliminated from the exec file. 

PMATRIXA and JMATRIXA contain the working "A" matrix 
for each of the programs. The working A matrix contains 
such information as the probability of success for each 
failure cause, the number of trials to failure for each 
cause, the system cause of failure, the phase number, and, 
in the case of JMATRIXA, the adjusted number of trials and 
adjusted number of failures. These FILEDEFS may be 
eliminated if desired. 

PREGMAT and JREGMAT contain the parameters computed in 
the REG, WREG1, WREG2, and WREG3 arrays as described on page 
2 of the program heading, (Appendix B, section 3). These 
matices pertain to the exponential regression and Weighed 
exponential regression models. These files are not 


required. 
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Ystar LISTING is produced only by DRG Fortran Al. Ee 
currently lists (arbitarily) the first 219 values for the 
tenth phase of testing. It is included merely as a 
verification of the implementation of this parameter. This 
FILEDEF may be eliminated. 

TRIALS DATA Al is also only produced by DRG Fortran Al. 
It lists the adjusted trials to failure for each phase. [In 
addition, the weights computed for use by methods 1 and 2 
are contained here. This FILEDEF is not required. 

EST OUT Al is produced only by JIMC Fortran Al and 
contains each phase estimate generated by the Woods Weighted 
Average Estimate Model for up to the first 5 simulations 
run. The remaining FILEDEFs with File Type "OUT" contain 
Similar information generated by the filename models. A 
file name beginning with the letter "J" corresponds to 
output produced by JIMC Fortran Al, while a file name 
starting with the letter "P" originated from DRG Fortran Al. 
These FILEDEFs should be retained. 

All of the remaining FILEDEFs correspond to files 
generated solely for plotting purposes. Each contains the 
mean reliability estimates of each phase in the Filetype 
"NUM" file and the corresponding standard deviation of these 
estimates in the Filetype "SDV" file. The information in 
these Files is duplicated in PRELIAB and JRELIAB LISTING. If 
the user is not interested in graphing results, these 24 


FILEDEFS may be eliminated. As explained earlier, the 
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integers in the filenames of these files correspond to the 
underlying actual reliability growth pattern. 
3. The Input Data File 

These instructions should be used in conjunction with a 
sample input file (See Appendix B, section 2) to aid in 
understanding. The easiest method for preparing an input 
file is to edit an existing input file. Due to the 
formatting of the "READ" statements in the simulation, it is 
imperative that all inputs be entered in the correct 
sequence. The simulation is built to read input from device 
number 10. 

The following steps should be allowed in order to 
produce an acceptable input file. All entries must be ona 
separate line in the input file although they may be 
anywhere in the line as long as the data entry is the first 
item encountered. 

1. Determine how many failure causes will be allowed in 
this simulation. This number must be an integer 
greater than or equal to one. There is no set limit 
on how large a number is possible. The capacity of 
the machine on which the program is being run will 
have some effect. Enter this number on the top line. 

2. Determine how many test phases are desired. Again, 
this number must be an integer greater than one. 
Enter this number on line 2. 

3. If the fixed phase reliability mode is desired than 


enter a loon line 3. If the constrained random growth 
mode is required than enter a O on line 3. 


ors) 


Enter the number of failures that will be allowed in 


Phase 1 on line 4. This number must also be an 
integer. A test-fix-test scenario would have one 
failure per phase, for example. Repeat this process 


on successive lines until you have entered the number 
of failures allowed for each phase that will be 
tested. Remember, if you defined the test as being 10 
phases, then you should have 10 separate entries, each 
on its own line, for this step. 


If the constrained random growth option is’ selected 
(you should have entered a O at line 3) then you must 
now enter the probability of success due to each 
failure cause for the first phase. This step should 
contain as many lines as the number of failure causes 
you identified in Step 1. The probability of success 
due to a failure cause is simply the probability of 
the failure cause occurring subtracted from one. 


If the fixed phase reliability option is selected (you 
should have entered a1 at line 3) then the process of 
entering the probabilities of success is a little 
more involved. Fixed phase reliability means that one 
is controlling the actual system reliability at each 
phase of the test. Since this system reliability is 
merely the product of the probabilities of success of 
all the failure causes at each phase then the user 
must enter these probabilities. These probabilities 
are entered by failure cause for each phase. If there 
are two failure causes and 3 phases in a test then the 
first entry will be the probability of success due to 
failure cause 1 in phase 1; the next entry will be the 
probability of success due to failure cause 1 in phase 
2 and then failure cause 1 in phase 3. After failure 
cause one has been entered for all three phase then 
failure cause 2's probabilities of success should be 
entered by by phase. Remember that each data entry 
must be on its own separate line. If X is the number 
of failure causes and Y is the number of test phase 
then this step should result in XY total data input 
lines. This number does not alter if more than one 
failure is permitted in a phase. 


re) 


LOe 


i les Le 


The next item that must be entered concerns use of the 
standard failure discount method. Even if you desire 
to use the C.L. discount method or if you do not 
desire to discount previous failures at all you must 
still place a value in this line. This line requires 
you to enter the discount interval or the number of 
successful trials that must occur between applications 
of the standard discount method. This number should 
be an integer. 


The next item required is the discount fraction. This 
value also applies to the standard discount method. 
If you do not want to discount previous failures then 
your should enter 0.0 on this line. If you desire to 
discount previous failures by 50 percent each time the 
discount method is applied then you should enter 0.50 
on this line. If you are using the C.L. discount 
method, you must still enter a value here, although it 
will not be used by the models in estimating 
reliability. 


The next required entry is the random number seed for 
random number ggnerator. Any number greater than zero 
and less than 2 will suffice. 


The FRIMP must be entered next. This value is only 
used if the constrained random growth mode is selected 
at line 3. This number must be between O and 1 and 
represents the fraction by which reliability will 
improve from phase to phase. It is applied to the 
probabilities of success of the failure causes each 
time they cause system failure in a phase. Thais 
method is intended to represent repairs or 
improvements in the system during the test. An entry 


of O will result in constant reliability while entries 
close to one will result in rapid reliability growth. 
Even if you are using the fixed phase reliability 
option you must still enter a value here although it 
will not play any role in the simulation. 


The next item that must be entered is the number of 
replications desired. For the purposes of the thesis 
associated with this paper 500 replications were done 
for each reliability growth pattern. 
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BPS 


14. 


aye 


6 . 


The next 7 (8, for use with JIMC Fortran Al) lines of 
required input deal with intermediate and specified 
Output. They are all binary (O or 1) option 
statements. The first addresses intermediate output. 
If intermediate phase by phase data is desired then 
the user should enter a 1 here. If intermediate 
output is not desired then enter a QO. This output 
will get voluminous rather quickly, particularly if a 
large number of relications is requested. The 
remaining options deal with the estimates generated 
by each model at each phase. The final output (which 
will be given regardless of the options chosen here) 
represents the average of all these estimates. If 
phase estimates are desired then a 1 should be 
entered; if not then enter a QO. The order of entry is 
weighted average model (for JIMC Fortran program 
only), MLE with failure discounting model, MLE single 
phase model, exponential regression model, and 
weighted exponential regression methods 1 through 3. 
Again, each data entry must be on a separate line. 


The next required entry is the failure discounting 
option. If you desire to use the standard discount 
method then enter al. If you desire not to discount 
at all enter ail and make sure you have entered 0.0 as 
instructed at Step 8 above. If you desire to use the 
C.L. method of failure discounting then enter a 2 in 
this space. 


Next enter the value of the C.L. method parameter (the 
confidence interval). This number should be between O 
and 1 and will typically be in the range of 0.8 to 
0.99. 


The next item that must be provided is the C.L. 
discount interval. Normally, the C.L. method does not 
use a discount interval (defined similar to the 
definition of a discount interval for the standard 
discount method) but in the majority of cases applying 
this method with some type of specified interval will 
lead to better results. These last two items must be 
entered even if one is using a alternate method of 
discounting. The other method will be applied but 
values are required due to the particular formatting 
of the simulation. 


The final entries that must be made are the weights to 
be accorded each phase of development. Enter any 
fraction between zero and one for each phase on a 
separate line. Remember the sum of these fractions 
must equal one. These entries are required regardless 
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of whether the user desires to apply weighting Method 

3 or not. Note: Giving equal weight to all phase is 

equivalent to not weighting and will produce results 
identical to the exponential regression model. 

If the input file is established consistent with the 

format outlined above and with the sample input file 


provided in Section 2 of Appendix B then there should be no 


problem in obtaining results. 
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APPENDIX B: FILES AND PROGRAMS 


SAMPLE EXEC FILE 
SAMPLE INPUT DATA FILE 
DRG FORTRAN PROGRAM 
JIMC FORTRAN PROGRAM 


SAMPLE OUTPUT FROM DRG FORTRAN Al 
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&TRACE OFF 


vevededoiededededekdedehdedetehhhehkchekicletedeidtctckhcndkichkkechticdddecckhick nk hkehkRR IR 


&FN = DRG 
&FN1 = RELIAB 


&TYPE Do you need to compile your program ? (Y) 


1. SAMPLE EXEC File 


&READ VAR &R_COMPILE 

SIF 6RSCONPICE NEY “GOTO -KUN 
= FORDVS. SEN 
&IF &RC EQ 0 &GOTO -RUN 


&IYPE Your program did net comp? le, 
&TYPE Do you wish to view the program LISTING file? (Y) 


&READ VAR &RSP1 
SIF €RSP] EQ % BROWSE GEN LISTINGss 


&TYPE Do you wish to XEDIT the program file? (Y) 
&READ VAR &RESP1 


Gir SRESP IONE YY -cbxit oi 
&COMMAND XEDIT &FN FORTRAN A 


&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP2 


CIP ORESP 27 srs eGOlOu-h 


&EXIT 1 
-RUN 

FILEDEF 
Pieper 
PILEDEF 
PCE Dir 
PGE DEY 
PIGEDEF 
PIG VEE 
F LEE DEP 
PILE 
PILEDEE 
FILGDEDEE 
FILEDEF 
PT LEDEE 
PPbEDES 
FILEDEF 
FILEDEE 
BE i313 
FILEDEF 
PilibE DE 
PIU LEDE: 
FILEDEF 
PIED GE 
FILEDEF 
JOR bg 
FIGEPEE 
FIELEVEP 
POLED EE 
FILEDEF 
i Eso 8 


10 
g2 


DISK 
DISK 


. JWRES2 
. JWRES3 
. MLEWD1 


. MLEWDS9 


INPUT1 DATA Al 
Al NUM Al 


. AY NUM Al 


JRELEAB LIS PING Al DRE ess 
PRELIAB LISTING Al (LRECL 133 
JTHESIS OUT Al 


) PORES ES OUT Ad 
. PMATRIXA LISTING (LRECL 133 
. JMATRIXA LISTING (LRECL 133 


PREGMAT DATA Al 


. JREGMAT DATA Al 
. YSTAR LISTING (LRECL 133 


TRIALS DATA Al 


ribo) OUl Ad 


JMLEWD OUT Al 


. PMLEWD OUT Al 
. JMLESP OUT Al 
(VEMGESP OUlL-Ai 


JREGEST OUT Al 


~ PREGES! O07 Al 
. PWRES1 


OUT 
PWRES2 
PWRES3 
JWRES1 


MLEWD1 


check for errors: 





FILEDEF 53 DISK MLEWD9 SDV Al 

PILEDEF 72 DISK REGS NUM 41 

POLEDEF 71 DISK KEG6 SDV Al 

FILEDEF 74 DISK REG16 NUM Al 

FILEDEF 73 DISK REG16 SDV Al 

FOELEDEF 77 DISK M1P1 NUM Al 

FEILEDEF 76 DISK M1Pi SDV Al 

MELEDEF 79 DISK MZP1 NUM Ail 

FILEDEF 78 DISK M2P1 SDV Al 

FILEDEF 92 DISK M3P1 NUM Al 

FILEDEF 91 DISK M3P1 SDV Al 

PYEEDEF 18 DISK MiP9 NUM Al 

FILEDEF 17 DISK MiP9 SDV Al 

FILEDEF 94 DISK M2P9 NUM Al 

PRE DEF 935 DISK M2P9 SDV Al 

PELE DEF 96 DISK M3P9 NUM Al 

PEEe DEF 95 DISK M3P9 SDV Al 

POLEDEF 62 DISK MLESP1 NUM Al 

PILEDEF 61 DISK MLESP1 SDV Al 

FILEDEF 64 DISK MLESP9 NUM Al 

FILEDEF 63 DISK MLESP9 SDV Al 

FILEDEF 06 TERMINAL 

LOAD &FN (START 

&IF &RC EQ O &SKIP 9 

fie PE Your program did not run correctly; check for errors. 
&TYPE Do you wish to XEDIT the program file? (Y) 
& READ VAR &RESP3 

ee &RESP3 NE Y &EXIT 2 

&COMMAND XEDIT & FN FORTRAN A 

&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP4 

SeckESP4 EQ Y &GOTO -H 

S&EXIT 2 

eee YOUR OUTPUT IS IN THE FILE &FN1 LISTING A 
&TYPE Do you wish to BROWSE your output? (Y) 

& READ VAR &RESP 

&IF &RESP EQ Y &COMMAND BROWSE &FNi LISTING A 
SxrE Print your output file? (CY) 

&READ VAR &RESP7 

&IF &RESP7 EQ Y &COMMAND PRINT &FN LISTING A 
-REDO 

&TYPE Do you wish to XEDIT the program file? (Y/N) 
&READ VAR &RESP5 

&IF &RESP5 EQ Y XEDIT &FN FORTRAN A 

&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP6 

&RESP56 = &CONCAT OF &RESPS &RESP6 

SeeeRESP56 EQ YY &GOTO -H 

&IF &RESP6 EQ Y &GOTO -RUN 

& EXIT 
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2. SAMPLE INPUT DATA FILE 


NUMBER OF FAILURE CAUSES 
NUMBER OF PHASES ( NPHASE ) 


Odo ocoocoocoocoocCO CO OCC COO OCOD DOO OC COCO CO OOCO COCO OO OC OF KF RP RRP RR RRR eeu 


FIXED RELIABILITY OPTION ( 1: YES oes 

NUMBER OF FAILURES IN PHASE 1 

NUMBER OF FAILURES IN PHASE 2 

NUMBER OF FAILURES IN PHASE 3 

NUMBER OF FAILURES IN PHASE 4 

NUMBER OF FAILURES IN PHASE 5 

NUMBER OF FAILURES IN PHASE 6 

NUMBER OF FAILURES IN PHASE 7 

NUMBER OF FAILURES IN PHASE 8 

NUMBER OF FAILURES IN PHASE 9 

NUMBER OF FAILURES IN PHASE 10 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASE 1 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASE 2 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASES 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASE 4 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASE 5 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASE 6 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASE 7 
PROB. OF SUCCESS FROM CAUSE 1) IN Piast 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASE 9 
PROB. OF SUCCESS FROM CAUSE 1 IN PHASE 1 
PROB. OF SUCCESS FROM CAUSE 2 IN PHASE 1 
PROB. OF SUCCESS FROM CAUSE 2 IN PHASE 2 
PROB. OF SUCCESS FROM CAUSE 2 IN PHASE 3 
PROB. OF SUCCESS FROM CAUSE 2 IN PHASE 4 
PROB. OF SUCCESS FROM CAUSES 2) INerisscee. 
PROB. OF SUCCESS FROM CAUSE 2 IN PHASE 6 
PROB. OF SUCCESS FROM CAUSE 2 IN Phaske 
PROB. OF SUCCESS FROM CAUSE 2 IN PHASE 8 
PROB. OF SUCCESS FROM CAUSE 2 IN PHASEe 
PROB. OF SUCCESS FROM CAUSE 2 IN iasie! 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASE 1 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASE 2 
PROB. OF SUCCESS FROM CAUSEU3 1) Sern astm. 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASE 4 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASE 5 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASE 6 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASE 7 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASE 8 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASES 
PROB. OF SUCCESS FROM CAUSE 3 IN PHASE 1 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 1 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 2 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 3 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 4 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 5 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 6 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 7 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 8 
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Ooroo ooo Qa qo ooo © 


oO oOo OO OC OC OF OF oOo oO 0 0 a 


PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 9 
PROB. OF SUCCESS FROM CAUSE 4 IN PHASE 10 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 1 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 2 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 3 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 4 
PROB] OReoUCCESS ERO CAUSE 5 IN PHASE 5 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 6 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 7 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 8 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 9 
PROB. OF SUCCESS FROM CAUSE 5 IN PHASE 10 
NUMBER OF TRIALS AFTER FAILURE BEFORE A DISCOUNT IS APPLIED 


FRACTION EACH FAILURE IS DISCOUNTED 
RANDOM NUMBER SEED FOR GGUBFS UNIFORM (0,1) GENERATOR 
FRACTION RELIABILITY IMPROVES AFTER FAILING IN A PHASE 
NUMBER OF DESIRED REPETITIONS FOR THE SIMULATION 
INTERMEDIATE INPUT OPTION(C1: 


SAVE ALL MLE W/ DISCOUNTING ESTIMATES (1: 


iio ONO Nie OUTPUT ) 


SAVE ALL MLE SINGLE PHASE ESTIMATES (1: 
SAVE ALL UNWT'D REGRESSION ESTIMATES (1: YES; 0: NO ) 


TES; 
YEO: O: 


O: 7 NOw) 
NO ) 


Saye lie iOuma Wi DeREGw ESI IMATES (1: YES; 0: NO ) 
SAVE ALL METHOD 2 WT'D REG. ESTIMATES (1: YES; 0: NO ) 
SAVE ALL METHOD 3 WT'D REG. ESTIMATES (1: YES; 0: NO ) 


DISCOUNTING OPTION (1: 


PERCENT C. I. 


C.L. DISCOUNT INTERVAL 


WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 
WEIGHT 


FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 


PHASE 


1 


PHASE 2 


PHASE 
PHASE 
PHASE 
PHASE 
PHASE 
PHASE 
PHASE 
PHASE 


h§ (© 00 NO F* NO 
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STRAIGHT % ; 


22 hoy Dp METHOR 
FOR C.L. DISCOUNTING METHOD(MUST HAVE A VALUE 


3. DRG FORTRAN PROGRAM 


Vevedesevesesedesedevedotedesesesedtovededesededesesesede Tedesco se fe hesevese de Reeves tedetevesedeheFedesetesesesetedesesedeFevesesete HARRI 


Fob th oh ok ok ok ot 


DISCRETE RELIABILITY GROWTH SIMULATION 
PROGRAMMED BY : JAMES E. DRAKE, 
JAMES D. CHANDLER, 
AND PAM A. MARKIEWICZ 


LAST MODIFIED 11 JUN 1988 


THE FOLLOWING EXTERNAL FILES ARE USED BY THE PROGRAM 


INPUT 


THESIS. 


RELIAB: 


MLEWD 


MLESP 


REGEST : 


WREST1 : 


WREST2 : 
WREST3 | 


DATA AND PARAMETER INPUT FILE (DEVICE # 10) 
OUTPUT FILE CONTAINING INTERMEDIATE COMPUTATIONS 
(DEVICE 725) 


OUTPUT FILE CONTAINING FINAL RESULTS OF THE SIMULATION 


(DEVICE # 35) 

OUTPUT FILE CONTAINING MLE ESTIMATES USING DISCOUNTING 
FOR EACH PHASE AND EACH REPLICATION 

(DEVICE # 55) 

OUTPUT FILE CONTAINING MLE ESTIMATE FOR EACH SINGLE PHASE 
AND ALL REPLICATIONS USING NO DISCOUNTING 

(DEVICE # 65) 

OUTPUT FILE CONTAINING EACH PHASE ESTIMATE FOR EACH 
REPLICATION OF THE EXPONENTIAL REGRESSION ESTIMATE 
(DEVICE 7 75) 

OUTPUT FILE CONTAINING EACH PHASE ESTIMATE FOR EACH 
REPLICATION OF THE WEIGHTED EXPONENTIAL REGRESSION 
ESTIMATE USING WEIGHTING METHOD 1) CDE tei) 

SAME AS ABOVE USING WEIGHTING METHOD 2 (DEVICE #39) 
SAME AS ABOVE USING WEIGHTING METHOD 3 (DEVICE #49) 


THE FOLLOWING IS A LIST OF KEY ARRAYS USED IN THE SIMULATION 


A 


NFAPH 


NFCAUS : 


NTRIAL : 


TRIALS : 


Pakeol = 


MAIN WORKING ARRAY CONTAINS PROBABILITY OF SUCCESS FOR 
EACH FAILURE CAUSE, NUMBER OF TRIALS UNTIL FAILURE FOR 
EACH FAILURE CAUSE AND THE SYSTEM, CAUSE OF FAILURE, 
PHASE NUMBER, ADJUSTED NUMBER OF TRIALS AND ADJUSTED 
NUMBER OF FAILURES 

DIMENSION ( ((2**#CAUSES)+6) ,#FAILURES ) 

CONTAINS THE NUMBER OF FAILURES IN EACH PHASE 
DIMENSION (1,#PHASES) 

BINARY ARRAY USED TO DETERMINE 
A PHASE 

DIMENSION ( 1,/#FAILURE CAUSES) 
CONTAINS THE NUMBER OF TRIALS SINCE LAST FAILURE OR 
DISCOUNTING FOR EACH FAILURE CAUSE 

DIMENSION ( 1,/FAILURE CAUSES ) 

CONTAINS THE TOTAL NUMBER OF ADJUSTED TRIALS IN A PHASE 
DIMENSION ( 1, # OF PHASES) 

RECORDS THE PHASE ESTIMATE FOR EACH ESTIMATOR WITHIN A 
SINGLE REPLICATION 


IF A FAILURE OCCURRED IN 
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a a a a a a a a 


« af. 
pS wv 


eo of ob tt 


REG 


WREG1 


WREG2 
WREG3 


DIMENSION (6,#PHASES ) 

ROW 1: MLE WITH DISCOUNTING 

ROW 2 : SINGLE PHASE MLE 

ROW 3 : WOODS REGRESSION ESTIMATE (UNWEIGHTED) 

ROW 4 : WEIGHTED REGRESSION ESTIMATE (METHOD 1) 

ROW 5 : WEIGHTED REGRESSION ESTIMATE (METHOD 2) 

ROW 6 : WEIGHTED REGRESSION ESTIMATE (METHOD 3) 
CONTAINS ACTUAL SYSTEM RELIABILITY IN EACH PHASE 
DIMENSION (1,#PHASES) 

CONTAINS YJK VALUES UP TO 1000 

DIMENSION (1,1000) 

CONTAINS THE YK VALUES UP TO 1000 

DIMENSION (#PHASES, 1000) 

CONTAINS THE VARIANCE OF THE ESTIMATE AT EACH PHASE 
DIMENSION (1,/#PHASES) 

CONTAINS THE RUNNING SUM OF 1/VAR(K) UP TO THE KTH PHASE 
DIMENSION (1,#PHASES) 

CONTAINS THE RUNNING SUM OF VAR(K) UP TO THE KTH PHASE 
DIMENSION (1,/#/PHASES) 

CONTAINS THE RUNNING SUM OF WEIGHTS UP TO THE KTH PHASE 
DIMENSION (1,/#PHASES) 


CONTAINS THE WEIGHT GIVEN TO EACH RELIABILITY ESTIMATE AT 


EACH PHASE. DIMENSION (3,#PHASES) 
ROW 1: WEIGHT = (1/VARCK))/(SUM OF 1/VAR(K) UP TO 
THE KTH PHASE) 
ROW 2: WEIGHT = VAR(K)/(SUM OF VARCK) UP TO KTH PHASE) 
ROW 3: WEIGHT = USER INPUTTED VALUES. 


ARRAY USED TO COMPUTE THE EXPONENTIAL REGRESSION ESTIMATE 


DIMENSION (5,#PHASES) 

ROW 1: K BAR 

ROW 2: Y BAR 

ROW 3: Y BAR FOR THE PHASE 
ROW 4 : B HAT 

ROW 5: A HAT 


ARRAY USED TO COMPUTE THE WEIGHTED EXPONENTIAL REGRESSION 


ESTIMATE USING METHOD 1 WEIGHTS. DIMENSION (4,##PHASES) 
ROW 1: X BAR(W) ...(A.K.A. K-BAR(W)) 

ROW 2: B HAT(W) 

ROW 3: A HAT(W) 

ROW 4: Y BAR(W) 


SAME AS ABOVE USING METHOD 2 WEIGHTS 
SAME AS ABOVE USING METHOD 3 WEIGHTS 


ele REMAINING ARRAYS ARE USED TO COMPUTE THE MEAN AND VARIANCE 
OF EACH ESTIMATE AT EACH PHASE. THEY ALL HAVE THE SAME DIMENSIONS 
AND STRUCTURE 


MLEWD 
MLESP 


REGE Ol ; 
WREST1 : 


DIMENSION (4,#PHASES ) 

ROW 1 : RUNNING SUM OF ESTIMATES 

ROW 2 : RUNNING SUM OF SQUARED ESTIMATES 
ROW Gon AN OF THE ESTIMATES 

ROW 4 : STANDARD DEVIATION OF THE ESTIMATES 


VALUES FOR THE MLE WITH DISCOUNTING 

VALUES FOR THE SINGLE PHASE MLE 

VALUES FOR THE EXPONENTIAL REGRESSION ESTIMATE 

VALUES FOR THE WEIGHTED EXPONENTIAL REGRESSION ESTIMATE 
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* USING METHOD 1 WEIGHTS. % 
* WREST2 : VALUES FOR THE WEIGHTED EXPONENTIAL REGRESS UGhee ee * 


ve USING METHOD 2 WEIGHTS. = 
* WREST3 : VALUES FOR THE WEIGHTED EXPONENTIAL REGRESSIONS Ss Pitan a 
¥ USING METHOD 3 WEIGHTS. - 


Cate nla elon 


C DEFINE AND DIMENSION VARIABLES 


PARAMETER (NR=50,NC=200) 

INTEGER REP,DISOPT,FRELOP,LDI,ALD 

REAL"4 MIN 

REAL*8 DSEED,MLESP ,MLEWD,EUL,SUM1,SUMS1,SUM2,SUMS2,SUM3 
C,SUMS3,SUM1A 

DIMENSION NFAPH(NR),A(NR,NC) ,NFCAUS(NR) ,NTRIAL( NR) , PHREST(6,NR), 
CMLEWD(4,NR) ,MLESP(4,NR) ,REGEST(4,NR) ,AREL(NR) , YIK(1000), 
CREG(5,NR) , YSTAR( 50,1000) , TADJTP( 1000) , VAR(NR) , SVAR1(NR) ,W(3,NR), 
CWREG1(4,NR) ,WREST1(4,NR) ,WREST2(4,NR) ,WREST3(4,NR), 
CSVAR2(NR) ,WREG2(4,NR) ,WREG3(4,NR) , SUMW(NR) 


C READ IN THE NUMBER OF CAUSES TO BE USED ( NCAUSE ) AND THE NUMBER 
C OF PHASES ( NPHASE ) IN THE TEST 


READ(10,**) NCAUSE 
READ(10,**) NPHASE 


C CHECK IF FIXED RELIABILITY OPTION IS DESIRED. FIX EULER'S NUMBER. 


READ(10,*) FRELOP 
EUL = 0.5772156648 


CREATE VARIABLES FOR THE ROW INDICES OF THE WORKING MATRIX ( A ) 
IPHASE: PHASE 
ISYSPR: ACTUAL COMPONENT RELIABILITY 
INTR: NUMBER OF TRIALS UP TO AND INCLUDING FAILURE 
IFAILC: CAUSE OF THE FAILURE 
IADJF: ADJUSTED NUMBER OF FAILURES ED 
AFTER DISCOUNTING HAS BEEN APPLIED 
IADJT: ADJUSTED NUMBER OF TRIALS AFTER DISCOUNTING HAS BEEN APPLIED 


Ge a gp a Sa Dl ep Fa J SP 


IPHASE = (2*NCAUSE)+1 
ISYSPR = IPHASE +1 
ENTR = 1Si¢SER e+ 1 
IFAILC = INTR + 1 
IADJF = IFAILC + 1 
IADJT = IADJF + 1 


READ IN THE NUMBER OF FAILURES IN EACH PHASE ( NFAPH(I) ) AND 
COMPUTE THE TOTAL NUMBER OF FAILURES IN THE TEST ( NFAIL ) 


So ee 


NFAIL = 0 
DO 10 I=1,NPHASE 
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C 
C 


C 


C 


ac) Ge) Gs CG) OO a aaa 


READ(10,*) NFAPH(T) 
NFAIL = NFAIL + NFAPH(1) 


10 CONTINUE 


1/7 


16 
Iles, 


20 


INPUT THE PROBABILITY OF SUCCESS IN A SINGLE TRIAL FOR EACH CAUSE 
IN EACH PHASE IF FRELOP EQUALS ONE. 


IF (FRELOP .EQ. 1) THEN 
DO 15 I=1,NCAUSE 
L=1 
DO 16 J=1,NPHASE 
READ(10,*) QQ 
DO 17 K= L,L+NFAPH(J)-1 


ACI, K) = QQ 
CONTINUE 
= L + NFAPH(J) 
CONTINUE 
CONTINUE 
ELSE 


INPUT THE PROBABILITY OF SUCCESS IN A SINGLE TRIAL FOR EACH CAUSE 
iN THE FIRST PHASE IF FRELOP EQUALS ZERO. 


DO 20 I=1,NCAUSE 
READ(10,**) A(I,1) 
CONTINUE 


ENDIF 


INPUT THE REMAINING VARIABLES , THE NUMBER OF SUCCESSFUL TRIALS 
BEFORE A DISCOUNT IS APPLIED (N); THE DISCOUNT FACTOR (R); THE SEED 
FOR THE RANDOM NUMBER GENERATOR, GGUBFS, (DSEED); RELIABILITY 
GROWTH FRACTION (FRIMP); TRIGGER FOR PRINTING INTERMEDIATE OUTPUT 
LOPT ) 
TRIGGERS FOR SAVING EACH ESTIMATE AT EACH PHASE FOR EACH ESTIMATOR 
TOPT1 : MLE WITH DISCOUNTING 
Ori 2). SINGLE PHASE MLE 
POEl oe se XPONENTIAL REGRESSION MODEL 
IOPT4 : WEIGHTED EXPONENTIAL REGRESSION MODEL (METHOD 1) 
TOPT5 : WEIGHTED EXPONENTIAL REGRESSION MODEL (METHOD 2) 
IOPT6 : WEIGHTED EXPONENTIAL REGRESSION MODEL (METHOD 3) 
DISCOUNTING OPTION TRIGGER (DISOPT); LLOYD FAILURE DISCOUNTING 
PARAMETER (GAMMA); LLOYD DISCOUNT INTERVAL 


READ(10,*) N 

READ(10,*) R 

READ(10,**) DSEED 
READ(10,*) FRIMP 
READ(10,*) NREP 
READ(10,*) IOPT 
READ(10,*) IOPT1 
READ(10,*) IOPT2 
READ(10,*) IOPT3 
READ(10,**) IOPT4 


wh 


READC 107%) 1@EiS 
READ(10,*) IOPT6 
READ(10,*) DISOPT 
READ(10,*) GAMA 
READ(10,*) LDI 


Cs=sec READ IN THE USER INPUTTED WEIGHTS FOR COMPUTATION OF METHOD: 
DO 25 I=1,NPHASE 
READ(10,*) W(3,1) 
25 CONTINUE 


XNREP = NREP 
DSEED1 = DSEED 


C INITIALIZE THE ARRAYS USED TO COMPUTE THE MEAN AND STANDARD DEVIATION 
C OF EACH ESTIMATOR 


DO 30 J=1,NPHASE 

DO 30 I=1,4 
MLEWD(I,J) 
MLESP(I,J) 
REGEST(1,J) 
WREST1(I,J) 
WREST2(I,J) 
WREST3(I,J) 


Oo © 
Oo. eo Oo = 
« s « « Cc & 


Sl DE eS 


30 CONTINUE 


DO 31 J=1,NPHASE 
DO- 31 = 6 


PERE SIT CL, J) 0.0 


oo CONTINUE 


YoK@l) — 2080 
DO 35 I=1,999 
YIK(I+1) = YJK(1I) + 1.0/1 
35. CONTINUE 


C COMPUTE AND STORE THE YSTAR VALUES UP TO 1000 
DO 39 J=1,NPHASE 
YSTAR(J,1) = 1. /NFAPH(J) 
DO 41 I = 1,999 
YSTAR(J,I+1)= YSTAR(J,I) + 1.0/(NFAPH(J)+1) 
41 CONTINUE : 
39 CONTINUE 


C PRINT OUT THE YSTAR MATRIX FOR THE FIRST 219 ENTRIES OF THE 1OTH PHaSE 


WRITE (90,*) (YSTAR(10,J), J=1,219) 


C COMPUTE AND STORE K BAR FOR THE EXPONENTIAL REGRESSION MODEL 


SUM = 0.0 
DO 50 I=1,NPHASE 


SUN eS Cee Ih 
REGO Ig) e= SUN it 
50 CONTINUE 


C MAJOR REPETITION OF THE SIMULATION LOOP 
DO 55 REP=1,NREP 


52 CONTINUE 
DO 500 REP=1,NREP 


C INITIALIZE FAILURE CAUSE VECTOR (NFCAUS) 
C COMPUTE THE INITIAL SYSTEM RELIABILITY 


REL = 1. 
DO 60 I=1,NCAUSE 
NFCAUS(I) = 0 
REL = REL * A(I,1) 
60 CONTINUE 


C INITIALIZE COLUMN (FAILURE # ) COUNTER FOR THE WORKING ARRAY (A) 
J=1 


C LOOP TO COMPUTE THE NUMBER OF TRIALS UP TO AND INCLUDING FAILURE 
C AND THE CAUSE OF FAILURE FOR EACH FAILURE IN EACH PHASE 


DO 130 K=1,NPHASE 


C SKIP ACTUAL COMPONENT RELIABILITY COMPUTATION AFTER FIRST REP 
C AND FOR FIRST FAILURE 


PCa. 1 )GolO 75 
PCR Giles GOlO.75 
nee 1. 


fee FIXED RELIABILITY OPTION IS SELECTED THEN PHASE RELIABILITIES 
C ARE COMPUTED AS FOLLOWS 


IF (FRELOP .EQ. 1) THEN 
DO 65 I=1,NCAUSE 
REL = REL*A(I,J) 
NFCAUS(I) = 0 
65 CONTINUE 
ELSE 


C COMPUTE NEW ACTUAL RELIABILITY FOR THE COMPONENT IN PHASE K 
DO 70 I=1,NCAUSE 


C INCREASE CAUSE PR(SUCCESS) IF IT CAUSED FAILURE IN THE PREVIOUS PHASE 
C COMPUTE NEXT PHACE RELIABILITY AND REINITIALIZE NFCAUS (NOT USED IF 
feet DXED PHASE RELIABILITY OPTION IS SELECTED). 


i(k 


PECNECAUS@I 2 EO. 1) THEN 

ACI,J) = ACI,(J-1)) + (C1. - ACL, (JI ee rER een 
EGoBIFCNECAUS( Tj} NE. 1) iia 

(led l= ACL, (J-ie 

Book 

ENDIF 

REL = REL*A(I,J) 

NFCAUS(I) = 0 


70 CONTINUE 
ENDIF 

iis: J1 = 1 
TROT =] 020 


CGC COMPUTE THE NUMBER OF TRIALS UP TO AND INCLUDING FAILURE AND THE 
C CAUSE OF FAILURE FOR EACH FAILURE IN THE PHASE 


DO 120 L=1,NFAPH(K) 
TFCREESCT Goro. 20 
TFC JIS BOs GOTOms> 
IP (GRELOr-. 20s) COnlenes 
DO 80 I=1,NCAUSE 

AGP = ACT (Joi) 


80 CONTINUE 
85 AC ISYSPR,J) = REL 

AC(IPHASE,J) = K 
90 MIN = 7. 2E75 


DO 110 I=1,NCAUSE 
C ASSIGN # TRIALS FOR CAUSES WITH PR(SUCCESS) = 0 OR 1 


IF(AC(I,J).GE.1.) THEN 
A((I+NCAUSE) ,J) = 7.2E75 
GOTO 100 

ELSEIF(A(I,J).EQ.0.) THEN 
ACCI+NCAUSE),J) = 1. 


GOTO 100 
ELSE 
ENDIF 
C CONVERT UNIFORM (0,1) RANDOM VARIABLE TO GEOMETRIC (i$ TRIALS UNTIL 
C FAILURE ) FOR EACH FAILURE CAUSE. RECORD THE MIN # TRIALS FOR THE 
C CAUSES AS THE SYSTEM # TRIALS UP TO AND INCLUDING FAILURE AND 
C RECORD THE FAILURE CAUSE 
A(CI+NCAUSE) ,J) = INT(1.+(LOG(GGUBFS(DSEED))/LOG(A(I,J)))) 
100 IF(A(( I+NCAUSE) ,J). LE. MIN) THEN 
MIN = A((I+NCAUSE),J) 
IMIN = I 
ELSE 
ENDIF 
110 CONTINUE 


ACIFAILC,J) = IMIN 
NFCAUSCIMIN) = 1 
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C COMPUTE THE TOTAL # OF TRIALS FOR THE MLE SINGLE PHASE ESTIMATE AND 
C INCREMENT FAILURE # COUNTERS 


ACINTR,J) = MIN 
TREO — Sik atet ee AC INIR, J) 
dosed + cl 
ewes 1 
120 CONTINUE 


C COMPUTE THE MLE ESTIMATE OF COMPONENT RELIABILITY FOR THIS PHASE AND 
C COMPUTE THE RUNNING SUM OF ESTIMATES AND THE SUM OF ESTIMATES SQUARED 
C FOR COMPUTATION OF THE MEAN AND STANDARD DEVIATION OF THE ESTIMATE 


Pokeomee sk )o= (ik tO! = NFAPHCK))/TRIOT 
PEE SE (ek) = MEESPC( I,K) + PHRESI(2,1K) 
MIR Sriezek) = MEESPC2,K) E CPHREST CZ, K)* 2) 


130 CONTINUE 


C INITIALIZE THE ADJUSTED NUMBER OF FAILURES TO 1 AND THE COUNT OF THE 
C NUMBER OF TRIALS SINCE FAILURE OR DISCOUNTING (NTRIALS(I) ) TO 0 
C IN PREPARATION FOR THE DISCOUNTING ROUTINE 


DO 140 J=1,NFAIL 
ACIADJF,J) = 1. 
140 CONTINUE 


DO 150 I=1,NCAUSE 
NTRIALCI) = 0 
150 CONTINUE 


C DISCOUNTING ROUTINE REVIEWS ALL PAST FAILURES AND CAUSES TO DATE 
C AND DETERMINES IF THE DISCOUNTING CONDITIONS HAVE BEEN MET. COMPUTES 
C THE ADJUSTED FAILURES, THE ADJUSTED # OF TRIALS AND YJK 
C INITIALIZE THE TOTAL ADJUSTED TRIALS IN A PHASE VECTOR, TADJTP; 
C INITIALIZE THE SUM OF THE VARIANCES,SVAR,THE VARIANCES, VAR 
C AND THE WEIGHTS,W FOR THE WEIGHTED EXPONENTIAL REGRESSION MODEL. 
DO 155 I = 1, NPHASE 
ATG) = © 
SVAR1(I) = 0 
SVAR2(I) = 0 
SUMW(I) = 0 
VAR(I) = 0 
von a0 
W(2,1) = 0 


155 CONTINUE 


SUMIA = 0 
J = 0 
LL = 1 
DO 300 K=1,NPHASE 
DO 200 L=1,NFAPH(K) 
Je=— soo 1 


S 


GC UPDATES THE NUMBER OF TRIALS SINCE FAILURE OR DISCOUNTING FOR EACH 
C FAILURE CAUSE 


ICAUSE = INT(A(IFAILC,J)+. 5) 
DO 160 I=1,NCAUSE 
IF( ICAUSE. EQ. I) THEN 
NTRIAL(I) = 0 
ELSEIF(ICAUSE. NE. 1) THEN 
NTRIAL(I) = NTRIAL(I) + INT(ACINTR,J)+. 5) 
ELSE 
ENDIF 
160 CONTINUE 


200 CONTINUE 
C CHOOSE DISCOUNTING METHOD TO BE USED 
[EC DISORR NES SGele ae 
C PERFORM LLOYD'S FAILURE DISCOUNTING METHOD 


DO 170 I=1,J 
I1 = INT(AC(IFAILC,I)+. 5) 
IF(NTRIAL(1I1).EQ.0) THEN 
ACIDE 1) = ee 
Gone 170 
ELSE 
ENDIF 


C THIS 18 THE MODIFIED LLOYD METHOD USING A DISCOUNT INTERVaE is 
G ORIGINAL DISCOUNT METHOD MAY BE EMPLOYED BY SETTING LDDSiGsG. 


ALD = INT(NTRIAL(11)/LDI) 
IF(ALD .EQ. 0) THEN 
ACIADJF,I) = 1.0 
GO TO 170 
ELSE 
ACIADJF,I) = 1.0 - ((1. -GAMA)**(1.0/ALD)) 
ENDIF 


170 CONTINUE 
GOTO 210 


C PERFORMS STRAIGHT PERCENT FAILURE DISCOUNTING AND 
C COMPUTES THE ADJUSTED # OF FAILURES 


180 DOO T=159 
I1 = INTCACIFAILC,1)+. 5) 
TPONTRIALCTI) SECO) tan. 
ACIADJF,1I) = 1. 
ELSBIPCNIRTAL( 11) .GEa Nps tien 
ACIADJF,I) = ACIADJF,1)*((1. -R)**CNTRIALC1I1)/N)) 
ELSE 
ENDIF 
190 CONTINUE 
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C ADJUSTS THE # TRIALS SINCE FAILURE OR DISCOUNTING FOR THOSE CAUSES 
C THAT HAVE MET OR SURPASSED THE DISCOUNTING THRESHOLD 
C FOR THE STRAIGHT PERCENT DISCOUNTING METHOD 


DO 205 I=1,NCAUSE 
IF(NTRIAL(I).GE.N) NTRIAL(I) = MODCNTRIAL(I),N) 
205 CONTINUE 
210 TADJT = 0.0 


C COMPUTES THE ADJUSTED # OF TRIALS FROM THE ADJUSTED # OF FAILURES 
C AND COMPUTES THE SUM OF THE ADJUSTED # OF TRIALS FOR ESTIMATE COMP. 


C IF ADJUSTED FAILURES ARE APPROACHING 0 THEN ADJUSTED TRIALS MUST 
fe DE PRE-SET. 


0) DF) 1S 


IF(ACIADJF,I) .LE. .0000001) THEN 
A(IADJF,I) = .0000001 
ENDIF 


AC IADJT,I) = ACINTR,I)/ACIADJF,I) 
Tap = TAPE -) ACTADIT Ay 
240 CONTINUE 


C COMPUTE THE ADJUSTED NUMBER OF TRIALS IN A PHASE 
DO 245 M = LL,(LL+NFAPH(K) - 1) 
TADJTP(K) = TADJTP(K) + AC IADJT,M) 
245 CONTINUE 
LL = LL + NFAPH(K) 
C COMPUTE THE MLE ESTIMATE OF PHASE RELIABILITY USING DISCOUNTING 
PHREST(1,K) = (TADJT - J)/TADJT 
COMPUTE Y-BAR AND Y-BAR FOR THE PHASE USING THE YSTAR MATRIX IF 


THE ADJUSTED NUMBER OF TRIALS IS LESS THAN OR EQUAL TO 1000 
AND USING EULER'S APPROXIMATION IF THE VALUE IS GREATER THAN 1000. 


Gre? ©) 


IF ( NINTCTADJTP(K)) . EQ. NFAPH(K) )THEN 
REGCS,K) = 0.0 

ELSE IF (NINT(CTADJTP(K) - NFAPH(K) ) .LE. 1000 ) THEN 
REG(3,K) = YSTARCK, NINTCTADJTP(K) -NFAPH(K) ) ) 

ELSE 

IX = NINTCTADJTP(K)-1 ) 

X = IX 

Q 


i 
S = X+2 

Ce COG CD (2) — P/Q) -C1/( GTS) ) 
REG(3,K) = U - YJK(NFAPH(K)~-1) 

END IF 


77 


SUM1IA = SUM1A + REG(3,K) 
REG(2,K) = SUM1A/K 


CG COMPUTE THE EXPONENTIAL REGRESSION ESTIMATE BEGINNING WITH B HAT 


SUM = 0.0 
SUMS = 0.0 
IF (K.EQ.1) GOTO 252 
DO 250 I = 1,K 
SUM = SUM + ((I-REG(1,K))*REG(3,1)) 
SUMS = SUMS + ((I-REG(1,K))**2) 
250 CONTINUE 
REG(4,K) = SUM/SUMS 


G >. COMPUTE A HAT 


REGCS,K) REG(2,K) 9 =e CREG( 47h) REGCIe ap 


GC COMPUTE AND STORE THE EXPONENTIAL REGRESSION ESTIMATE 


PHREST(3,K) = 1.0 - EXP(-(REG(5,K) + (REG(4,K)*K))) 
IF(PHREST(3,K). LT. 0.0) PHREST(3,K)=0. 0 
GOTO 255 ~ 
257 PHREST(3,K) = 1.0 - EXP(-REG(3,1)) 
IF( PHREST(3,K). LT. 0.0) PHREST(3,K)=0. 0 


STORE THE RUNNING SUM OF THE ESTIMATES FOR THE CURRENT PHASE AND THE 
RUNNING SUM OF THE ESTIMATES SQUARED FOR COMPUTATION OF THE MEAN AND 
STANDARD DEVIATION OF EACH ESTIMATE FOR EACH RELIABILITY GROWTH 
MODEL 


C2) C1) 


255 MLEWD(1,K) = MLEWD(1,K) + PHREST(2,K) 
MLEWD( 2 ,K) MLEWD(2,K) + (PHREST(2,K)**2) 
REGEST(1,K) = REGEST(1,K) + PHREST(3,K) 
REGEST(2,K) = REGEST(2,K) + (PHREST(3,K)**2) 


C COMPUTE THE VARIANCES OF THE UNWEIGHTED EXPONENTIAL REGRESSION 
C ESTIMATES AND STORE THE RUNNING SUM OF THE VARIANCES FOR USE IN THE 
C WEIGHTED REGRESSION MODEL. 

VAR(K) = (NFAPH(K)**PHREST(1,K))/((1. -PHREST(1,K))**2) 

IF ( VAR(K) .LT. .0000001)THEN 

VAR(K) = .0000001 

END IF 

IF (XK GEOn Ge Tease 

SVAR1(K) = SVAR1(K-1) + (1. /VAR(K)) 

SVAR2(K) = SVAR2(K-1) + VAR(K) 

SUMW(K) = SUMW(K-1) + W(3,K) 

GO TO 259 

258 SVAR1(K) = (1. /VAR(1)) 
SVAR2(K) = VAR(1) 
SUMW(K) = W(3,1) 


78 


C COMPUTE THE WEIGHTS FOR EACH PHASE 


259 WC1,K) = (1. /VARCK))/SVAR1(K) 


261 


WC2,K) = VARCK)/SVAR2(K) 
C COMPUTE AND STORE X-BAR(W) AND Y-BAR(W) 


SUMX1 
SUMY1 
SUMA2 
SUMY2 
SUMX3 
SUMY3 
Doe Zou. 


| 
oCcoeo ee 


SUMX1 
SUMY1 
SUMX2 
SUMY2 
SUMXK3 
SUMY3 


CONTINUE 

WREG1(1,K) 
WREG1(4,K) 
WREG2(1,K) 
WREG2(4,K) 
WREG3(1,K) 
WREG3(4,K) 


1,K 


SUMXI1 
SUMY1 
SUMX2 
SUMY2 
SUMX3 
SUMY3 


(1. /VAR(K) )*1 

(1. /VAR(K) )*REG(3,I) 
VAR(K)*I 

VAR( K)*REG(3,I) 
W(3,K)*1 
W(3,K)*REG(3,1) 


SUMX1/SVAR1(K) 
SUMY1/SVARI1(K) 
SUMX2/SVAR2(K) 
SUMY2/SVAR2(K) 
SUMK3 / SUMW(K) 
SUMY3 / SUMW(K) 


C COMPUTE AND STORE B-HAT(W) 


AD) Me 


SUM1 
SUMS 1 
SUM2 
SUMS 2 
SUM3 
SUMS3 
JR CB INOR 


0 
0 
Oo: 
0. 
0 
0 
Q 


WO 2 7 el 
SUM1 
SUMS 1 
SUM2 
SUMS2 
SUM3 
SUMS3 

CONTINUE 


WREG1(2,K) 
WREG2(2,K) 
WREG3(2,K) 


Oo Ceo eo Oo 


iy Gonna) 272 


1, 


K 


SUM1 + (WC1,K)*( I1-WREG1(1,K))*REG(3,I)) 
SUMS1 + ( W(1,K)*((I-WREG1(1,K) )**2)) 
SUM2 + (W(2,K)**( 1-WREG2(1,K))**REG(3,1)) 
SUMS2 + ( W(2,K)*((1-WREG2(1,K) )**2)) 
SUM3 + (W(3,K)*( 1-WREG3(1,K))*REG(3,I)) 
SUMS3 + ( W(3,K)*((I-WREG3(1,K))**2)) 


SUM1/SUMS1 
SUM2/SUMS2 
SUM3/SUMS3 


C COMPUTE AND STORE A-HAT(W) 


WREG1(3,K) 
WREG2(3,K) 


WREG1(4,K) - C(WREG1(2,K)*WREG1(1,K)) 
WREG2(4,K) - (WREG2(2,K)*WREG2(1,K)) 


(e 


WREG3(3,K) = WREG3(4,K) - (WREG3(2,K)*WREG3(1,K)) 
C COMPUTE AND STORE THE WEIGHTED EXPONENTIAL REGRESSION ESTIMATE 


PHREST(4,K) 
PURE SI ork) 


1.0 - EXP( -(WREGI(3,K)?( WREGIC 2) ee 
1.0 - EXPC=(WREG2Z(3, K)+(WREG2( Zoho) tae) 


PHREST(6,K) = 1.0 - EXP(-(WREG3(3,K)+(WREG3(2,K)*K))) 
IF (PHREST(4,K) .LT. 0.0) PHREST(4,K) = 0.0 
IF (PHREST(5,K) .LT. 0.0) PHREST(5,K) = 0.0 
IF (PHREST(6,K) .LT. 0.0) PHREST(6,K) = 0.0 
GO TO 275 

272 PHREST(4,K) = 1.0 - EXP(-WREG1(4,1)) 
IF (PHREST(4,K) .LT. 0.0) PHREST(4,K) = 0.0 
PHREST(5,K) = 1.0 - EXP(-WREG2(4,1)) 
IF (PHREST(5,K) .LT. 0.0) PHREST(5,K) = 0.0 
PHREST(6,K) = 1.0 - EXP(-WREG3(4,1)) 
IF (PHREST(6,K) .LT. 0.0) PHREST(6,K) = 0.0 


C STORE THE RUNNING SUM OF THE WEIGHTED ESTIMATES FOR THE Cun aa® 

C PHASE AND THE RUNNING SUM OF THE ESTIMATES SQUARED FOR COMPUTATION OF 
C THE MEAN AND STANDARD DEVIATION OF EACH WEIGHTED ESTIMATE FOR THE 

C WEIGHTED EXPONENTIAL REGRESSION MODEL. 


275 WRESTI(1,K) = WRESTIC1,K) + PHREST Conk 
WRESTI(2,K) = WREST1(2,K)°+ CPRRES hee 2) 


WREST2(1,K) 
WREST2( 2 ,K) 


WREST2(1,K) + PHREST(5,K) 
WREST2(2,K) + (PHREST(5 ,K)**2) 


WRESI2¢1,4) 
WREOIGG2).K) 


WREST3(1,K) + PHREST(6,K) 
WREST3(2,K) + (PHREST(6,K)**2) 


C STORE THE ACTUAL PHASE RELIABILITY 
AREL(K) = ACISYSPR,J) 


C PRINT INTERMEDIATE OUTPUT IF REQUESTED AND THE NUMBER OF REPE DDR 
C IS NOT GREATER THAN 5 


PEC IOPT NEw) Gores 00 
IFCREF. Ol. Ss )=CUdOn 00 


WRITE(25,1000* REP,K 

1000 FORMAT(T16,'REPETITION NUMBER: ',14,' PHASE NUMBER: ‘ ,I4) 
WRITE( 25,1010) ACISYSPR,J) 

1010 FORMAT(22X,'ACTUAL COMPONENT RELIABILITY: ',F7.5) 
WRITEC 25, 1022). PHRESROlek) 

1022 FORMAT( 20X,'MLE ESTIMATE USING DISCOUNTING: ',F7.5) 
WRITE(25,1025) PHREST(2,K) 

1025 FORMAT(18X,'MLE ESTIMATE OF PHASE RELIABILITY: ',F7.5) 
WRITE( 25,1027) PHREST(4,K) 

1027 FORMAT(14X, ‘WEIGHTED REG. ESTIMATE ( METHOD 1 ) 5 Ee 5) 


SO 


Wk I TE@255 026) seHRESt(5 ,K) 


1026 FORMAT(14X,'WEIGHTED REG. ESTIMATE ( METHOD 2)  : ',F7.5) 
WRITE( 25,1029) PHREST(6,K) 
1029 FORMAT(14X,'WEIGHTED REG. ESTIMATE ( METHOD 3) =: ',F7.5) 


WRITE(25,1028) PHREST(3,K) 
1028 FORMAT(14X,'REGRESSION ESTIMATE OF PHASE RELIABILITY: ',F7.5) 
WRITE( 25,1030) 
1030 FORMAT(' ',' =) 
DO 260 I=1,NCAUSE 
WRITE(25,1035)1,ACI 
1035 FORMAT(12X, ‘CAUSE: 
CF10. 0) 
260 CONTINUE 
WRITE( 25 , 1036) 
1036 FORMAT(' ',' 
WRITE( 25,1040) 
1040 FORMAT(4X,'FAIL #',3X,'FAIL CAUSE' ,3X,'# TRIALS' ,3X,'ADJ # FAIL’ ,3 
CX,'ADJ # TRIALS') 
DOm20 lio 
MRE G25, 10501 AC lPAlucC. 1) ,ACINIR, I) AC LADJF,1),ACIADJT,1) 
MS OP PORMATIC 4X ,13,6X,F3.0,7X,F8.0,4X,F8.6,4X%,F12.0,3X,F11. 4) 
270 CONTINUE 
WRITE( 25 , 1060) 
1060 FORMAT(' ',///) 


J) ,ACCI+NCAUSE) , J) 


'13,'  PRCSUCCESS): ',F7.6,'  # TRIALS: ', 


300 CONTINUE 


Seen INTSEACH OF THE ESTIMATES TO THEIR APPROPRIATE OUTPUT FILE 
f IF REQUESTED 


401 IFCIOPT1.NE.1) GOTO 402 

WRITE(55,2000) (PHREST(1,1), I=1,NPHASE) 
402 IF(IOPT2.NE.1) GOTO 403 

WRITE(65,2000) (PHREST(2,1), I=1,NPHASE) 
403 IF(IOPT3.NE.1) GOTO 404 

WRITE( 75,2000) (PHREST(3,I), I=1,NPHASE) 
404 IF(IOPT4.NE.1) GOTO 405 

WRITE(15,2000) (PHREST(4,1), I=1,NPHASE) 
405 IF(CIOPT5.NE.1) GOTO 406 

WRITE( 39,2000) (PHREST(5,1), I=1,NPHASE) 
406 IF(IOPT4.NE.1) GOTO 500 

WRITE( 49,2000) (PHREST(6,1), I=1,NPHASE) 
2000 FORMAT(' ',30(F7. 6: 1X)) 


500 CONTINUE 


C PRINT OUT THE WORKING "A'’ MATRIX IN MATRIX LISTING Al 
DO 4050 J = 1,(2*NCAUSE)+6 
WRITE(81,*) (ACJ,1),I = 1,NFAIL) 

4050 CONTINUE 


C PRINT OUT THE THE TOTAL NUMBER OF TRIALS TO SYSTEM FAILURE IN EACH 
C PHASE 
DO 4051 J = 1,NPHASE 
WRITE(89,*) 'T ',TADJTP(J) 
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4051 CONTINUE 


C PRINT OUT THE WEIGHTS COMPUTED BY METHOD 1 TO THE TRIALS FILE 
DO 4052 J = 1,NPHASE 
WRITE(89,*) ‘M1 ',W(1,J) 
4052 CONTINUE 


C PRINT OUT THE WEIGHTS COMPUTED BY METHOD 2 TO THE TRIALS FILE 
DO 4053 J = 1,NPHASE 
WRITE(89,*) ‘M2 ',W(2,J) 
4053 CONTINUE 


C PRINT OUT THE REGRESSION PARAMETERS IN MATRIX FORM TO THE REGMAT FILE 


DO 4054 J=1,5 
WRITE(87,*) (REG(J,1),1 = 1,NPHASE) 
4054 CONTINUE 


C PRINT OUT THE WEIGHTED REGRESSION PARAMETERS IN MATRIX FORM TO THE 
C REGMAT FILE 
DO 4055 J = 1,4 
WRITE(87,*) (WREG1(J,1I),I 
4055 CONTINUE 


1,NPHASE) 


DO 4056 J = 1,4 
WRITE(67,**) (WREG2Z0), 1), 12= NPT ASE. 
4056 CONTINUE 


DO 4057 J = 1,4 
WRITE(87,*) (WREG3(J,1),I = 1,NPHASE) 
4057 CONTINUE 7 


C UPON COMPLETION OF ALL REPETITIONS, COMPUTE THE MEAN AND STANDARD 
C DEVIATION OF EACH ESTIMATE FOR EACH PHASE SKIPPING COMPUTATIONS IF 
C ONLY ONE REPETITION IS REQUIRED 


TE OONKEP, ihe) -GOnorool 
DO 600 I=1,NPHASE 


MLEWD(3,I) = MLEWD(1,1)/XNREP 
MLESP(3,I) = MLESP(1,1)/XNREP 


REGESI(3,1) = REGESI( 1.1) Nee 
WREST1(3,1) = WRESTIC1(1)/ ANREP 
WREST2(3,1) = WREST2(1,1)/XNREP 
WREST3(3,1) = WREST3(1,1)/XNREP 
MLEWD(4,1) = SQRT( (MLEWD( 2,1) -C(XNREP*(MLEWD(3,1)**2)))/(XNREP-1) ) 
MLESP(4,1) = SQRT((MLESP(2,1)-(XNREP*(MLESP(3,1)**2) ))/( XNREE Se 
REGEST( 4, I)=SQRT( CREGEST( 2, I )-(XNREP*“( REGEST( 3, 1)***2) ) )/ (ANREP Se 
WREST1(4, I)=SQRT( (WREST1(2, 1) -(XNREP*(WRESTIC3,1)**2) ))/( SNRED Se 
WREST2(4, 1)=SQRT( (WREST2( 2,1) -(XNREP*(WREST2(3,1)**2)))/(XNREP=1)) 
WREST3(4, 1)=SQRT( (WREST3( 2,1) -(XNREP*(WREST3(3, 1)**2) )) / CANRED ee 


600 CONTINUE 
C PRINT THE FINAL OUTPUT TABLE [O74 tie 


601 WRETECSS 23000) 


3000 FORMAT('0' ,T47,'DISCRETE RELIABILITY GROWTH SIMULATION’ ) 

WRITE( 35,3010) 
3010 FORMAT('-',T54,'MODEL PARAMETER SUMMARY’ ) 

WRITE( 35,3020) NCAUSE 
3020 FORMAT('0',T47,'NUMBER OF POSSIBLE FAILURE CAUSES ',14) 

IF (FRELOP .EQ. 1) GOTO 4000 

WRITE(35 , 3030) 
3030 FORMAT('0' ,T38,'CAUSE NUMBER' ,T64,'SINGLE TRIAL PR( SUCCESS ) FOR 

CPHASE 1') 

DO 3050 M=1,NCAUSE 

WRITE( 35,3040) M,AC(M,1) 
3040 FORMAT(' ',T43,12,T79,F8. 6) 
3050 CONTINUE 

WRITE( 35,3060) FRIMP 
3060 FORMAT('0' ,T37,'FRACTION CAUSE RELIABILITY IMPROVES AFTER FAILURE 

GSES. 6) 
5000 WRITE(35,3080) NPHASE 
3080 FORMAT('-',T48,'NUMBER OF PHASES IN THE SIMULATION ',I2) 

WRITE(35 , 3090) 
3090 FORMAT('0',T42,'PHASE NUMBER' ,T59,'NUMBER OF FAILURES IN THE FIRST 

C PHASE' ) 

DO 3110 M=1,NPHASE 

WRITE( 35,3100) M,NFAPH(M) 
SOO FORMATC (043 ,12,173,12) 
3110 CONTINUE 

WEeIME@S5 3120) NFAIL 
3120 FORMAT('0',T51,' TOTAL NUMBER OF FAILURES ',1I4) 

IF(DISOPT. EQ. 2) GO TO 3160 


LAN OSS 61.50) 
3130 FORMAT('-',T38,' DISCOUNTING PERFORMED USING THE CONSTANT FRACTION 
CMETHOD' ) 


WRITE( 35,3140) R 

3140 FORMAT('0' ,T44,' FRACTION EACH FAILURE IS DISCOUNTED ',F8. 6) 
WRITE( 35,3150) N 

3150 FORMAT(' ',T33,'NUMBER OF TRIALS AFTER A FAILURE BEFORE A DISCOUNT 
CALS WAPPLIBD 4) 
GO TO 3190 

3160 WRITE(35,3170) 

3170 FORMAT('-' ,T44,'DISCOUNTING PERFORMED USING THE LLOYD METHOD' ) 
WRITE( 35,3180) GAMA 

3180 FORMAT('0' ,T39,'PERCENT C.I. ( USED AS DISCOUNT FRACTION ) ',F8.6 
C) 
WRITE( 35,3185) LDI 

3185 FORMAT('0',T50,' LLOYD DISCOUNT INTERVAL: ',I3) 

3190 WRITE(35,3200) DSEED1 

3200 FORMAT('-',T46,'RANDOM NUMBER SEED USED ',F15. 2) 
WRITE(35,3210) NREP 

3210 FORMAT( ‘0° ,T37, 'NUMBER OF REPETITIONS OF THE SIMULATION PERFORMED 
Cel) 
WRITE(35 , 3220) 

3220 FORMAT('1' ,T61,'ESTIMATOR: ') 
WRITE( 35 , 3230) 

3230 FORMAT('0',T48,'SINGLE PHASE MLE WITHOUT DISCOUNTING' ) 
WRITE( 35 ,3240) 

3240 FORMAT('-' ,T60,'MEAN' ,T83, 'ESTIMATE' ,T109,'95 %') 
WRITE(35 , 3250) 
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3250 FORMAT(' ',T12,'PHASE NUMBER’ ,T29,'ACTUAL RELIABILITY' ,T52, ‘PREDIC 
CTED RELIABILITY’ ,T78,'STANDARD DEVIATION’ ,T101,'CONFIDENCE INTERVA 
CE) 


GC COMPUTE G.i1. FOR SINGLE PHASESMEE 


DO 3270 M=1,NFHASE 
CI = (1. 96*MLESP(4,M))/SQRT( XNREP) 
CIU = MLESP(3,M) + CI 
CIL = MEESP(3.é) = el 
WRITE( 35,3260) M,AREL(M) ,MLESP(3,M),MLESP(4,M),CIL,CIU 
3260 FORMAT('0' ,T17,12,T34,F8.6,1T58,F8. 6,T82,F9. 6,199, '( " {8.6 
Gono 
3270 CONTINUE 
WRITE( 35 , 3220) 
WRITE(35, 3280) 
3280 FORMAT('0',T42,'MAX LIKELIHOOD ESTIMATE USING DISCOUNTED FAILURES' 
C) 
WRITE( 35,3240) 
WRITE(35 , 3250) 


C COMPUTE C.I. FOR MLE WITH DISCOUNTING 


DO 3290 M=1,NPHASE 

CI = (1. 96*MLEWD(4,M))/SQRT(XNREP) 

CIU = MLEWD(3,M) + CI 

CIL = MLEWD(3,M) - CI 

WRITE( 35,3260) M,AREL(M) ,MLEWD(3,M) ,MLEWD(4,M) ,CIL,CIU 
3290 CONTINUE 


WRITE( 35 , 3220) 
WRITE( 35 , 3320) 
3320 FORMAT('0',T43, REGRESSION ESTIMATE USING DISCOUNTED FAILURES') 
WRITE( 35 ,3240) 
WRITE( 35,3250) 


C COMPUTE C.I. FOR EXPONENTIAL REGRESSION ESTIMATES 


DO 3330 M=1,NPHASE 

CI = (1. 96*REGEST(4,M))/SQRT(XNREP) 

CIU = REGEST(3,M) + CI 

CIL = REGEST(3,M) - CI 

WRITE( 35,3260) M,AREL(M) ,REGEST(3,M),REGEST(4,M) ,CIL,CIU 
3330 CONTINUE 


Csesc'¢ WEIGHTED REGRESSION (METHOD 1) 
WRETE Coo. o2 20) 
WEE Gson 0 32h) 
3321 FORMAT('0' ,T43, ‘WEIGHTED REGRESSSION ESTIMATE (METHOD 1 )') 


WRITE( 35,3240) 
WRITE( 35,3250) 


C COMPUTE C.I. FOR WEIGHTED EXPONENTIAL REGRESSION ESTIMATES 


DO 3331 M=1,NPHASE 
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= (1. 96*WREST1(4 ,M) )/SQRTC XNREP ) 
ClUe=sWReo! 135, i) el 
Cil, = WREST1C3,¢M) - CI 
WEE soe o2 CO eM annem yeWRESTI03,M),WREST1(4,M),CIL,CIU 
3331 CONTINUE 


0 =66ee WEIGHTED REGRESSION (METHOD 2) 
WRITE( 35,3220) 
WietMin@s> 53522) 
3322 FORMAT('0',T43,'WEIGHTED REGRESSION ESTIMATE (METHOD 2) ') 


WRITE( 35, 3240) 
WRITE( 35, 3250) 


C COMPUTE C.I. FOR WEIGHTED EXPONENTIAL REGRESSION ESTIMATES 


DO 3332 M=1,NPHASE 
= (1. 96*WREST2(4,M))/SQRT( XNREP) 
CIU = WREST2(3,M) + CI 
CIL = WREST2(3,M) - CI 
WRITE( 35,3260) M,AREL(M) ,WREST2(3,M),WREST2(4,M) ,CIL,CIU 
3332 CONTINUE 


a WEIGHTED REGRESSION (METHOD 3) 
WRITE( 35, 3220) 
WRITE( 35 , 3323) 
3323 FORMAT('0',T43,'WEIGHTED REGRESSION ESTIMATE (METHOD 3) ') 


WRITE( 35 ,3240) 
WRITE( 35,3250) 


G COMPUTE C.1I. FOR WEIGHTED EXPONENTIAL REGRESSION ESTIMATES 


DO 3333 M=1,NPHASE 

CI = (1. 96*WREST3(4,M))/SQRT( XNREP ) 

Clu = WREST3(3,M) + CI 

Tite wkhohs( 3, iM) =. Cl 

WERE TE@25, 3200) M,ARELCM) WREST3(3,M),WRES1T3(4,M),CIL,CIU 
e259 CONTINUE 


WRITE( 35 , 3340) 
3340 FORMAT('1' ,T59, 'RECAPITULATION' //) 
WRITE( 35, 13350) | 
3350 FORMATC( ' ,T3, PHASE! Die cCnUsl -'2s. MEAN (138, EST Roe 
C'MEAN' 163° ‘EST! aloe MEAN' ooweEse ,1103, NEAN ,2l13; EST ) 
WRITE(35, 3360) 
3360 FORMAT( ' ,T11,' RELIAB’ ,T28, oe ,T38,' STD! PT Somes | 065s 50) 17 
C7,'PHASE' ,T88, STD' ,T103, "REG! wigs Sip ) 
WRITE(35, 3370) 
3370 FORMAT( | Si CreeRG. 4135, DEVIATION’ ,T53, ‘W/D! ,T60, DEVIATION’ 
S176, MLE’ Te5. ' DEVIATION’ PLLOe, usd 4b) 10, ‘DEVIATION’ ) 
WRITE(35, 3375) 
275 FORMAT( ,1T28, EST'/) 
DO 650 I=1,NPHASE 
WRITE(35 , 3380) I,AREL(1) ,WREST1(3,1) ,WREST1(4,1) ,MLEWD( 3,1) 
C,MLEWD(4,1I) ,MLESP(3,I) ,MLESP(4,1),REGEST(3,1),REGEST(4,I) 


3380 FORMAT('0' ,T4,13,T11,F7.6,926,F7.6,136,F/7. 6.15150 708 loaner 
C,1T76,F7.6,T86,F7.6,T101.F7- omer tence) 
650 CONTINUE 


Ca Sas PAGE 2 OF RECAPITULATION: 


WRITE( 35,3381) 
3381 FORMAT('1',T59,'RECAPITULATION CONT. '//) 


WRITE(35 , 3382) 

3382 FORMAT( '-" ,T3, ‘PHASE’ ,T11, ACTUAL’ ,T28, ‘MEAN’ ,T38, 'EST' T53 
C'MEAN' ,T63, EST',T78, MEAN’ ,T88, EST’ ,T103, MBEAN’ , Tio oo 
WRITE( 35,3383) 

3383 FORMAT(' ',T11, ‘RELIAB' T28, REG’ ,T38,'STD’ ,T53, METHOD = e-8 
C'STD' ,T77, METHOD’ ,T88, STD’ ,T103, METHOD (113 oop 
WRITE(35 , 3384) 

3384 FORMAT(' ',T28,'EST' ,T35,'DEVIATION' ,T53,‘'ONE' ,T60, DEVIATION’ 
C,T78, TWO’ ,T85, DEVIATION ,T103, THREE , T1110, DEVIATIC Ue 
DO 651 J=1,NPHASE 
WRITE( 35,3385) I,AREL(1),REGEST(3,1),REGEST(4,1),WREST1(3,1) 
C ,WREST1(4,1) ,WREST2(3,1) ,WREST2(4,1) ,WREST3(3,1) ,WREST3(4,1) 

3385 FORMAT('0' ,T4,13,T11,F7. 6,1T26,F7.6,1T36,F/.6,151,F7.6,06), Fee 
G,.176,F7.6,186,F7.6,1401,F/.6, biter) 


' 


G PRINT RELIABILITY ESTIMATES TO FILES FOR PLOTTING FURPGsE: 


WRITE(51,3400) MLEWD(4,I1) 
WRITE(52,3400) MLEWD(3,I1) 
WRITE(61,3400) MLESP(4,1) 
WRITE( 62,3400) MLESP(3,I1) 
WRITE(71,3400) REGEST(4,I1) 
WRITE( 72,3400) REGEST(3,1) 
WRITE( 76,3400) WREST1(4,1) 
WRITE( 77,3400) WREST1(3,1) 
WRITE( 78,3400) WREST2(4,1) 
WRITE( 79,3400) WREST2(3,1) 
WRITE(91,3400) WREST3(4,1) 
WRITE(92,3400) WREST3(3,1) 
WRITE( 82,3400) AREL(I1) 
3400 FORMAT('' ,F7. 6) 


651 CONTINUE 


GO TO 6000 


4000 WRITE(35 ,4010) 
4010 FORMAT(1X,//,T50,' FIXED PHASE RELIABILITY OPTION’ ) 
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WRITE(35 ,4020) 
4020 FORMAT('-' ,T38,'PHASE NUMBER’ ,T78,' ACTUAL RELIABILITY’) 
DO 4030 M=1,NPHASE 
WRITE(35,4040) M,AREL(M) 
4040 FORMAT('0' ,T41,12,T83,F8. 6) 
4030 CONTINUE 
GO TO 5000 
6000 CONTINUE 


SOP 
END 
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4. JIMC FORTRAN PROGRAMNI 


DISCRETE RELIABILITY GROWTH SIMULATION 


PROGRAMMED BY JAMES E DRAKE, 
JAMES D CHANDLER, 

AND PAM A MARKIEWICZ 

LAST MODIFIED 11 JUN 1988 


THE FOLLOWING EXTERNAL FILES ARE USED BY THE PROGRAM 
INPUT : DATA AND PARAMETER INPUT FILE (DEVICE # 10) 
THESIS : OUTPUT FILE CONTAINING INTERMEDIATE COMPUTATIONS 
(DEVICE # 20) 
RELIAB: OUTPUT FILE CONTAINING FINAL RESULTS OF THE SIMULATION 
CUEVIECh a 30) 
EST : OUTPUT FILE CONTAINING EACH PHASE ESTIMATE FOR EACH 
REPLICATION OF THE WOODS WEIGHTED AVERAGE ESTIMATE 
(DEVICE #40) 
MLEWD : OUTPUT FILE CONTAINING MLE ESTIMATES USING DISCOUNTING 
FOR EACH PHASE AND EACH REPLICATION 
(DEVICE 7 50) 
MLESP : OUTPUT FILE CONTAINING MLE ESTIMATE FOR EACH SINGLE PHASG@ 
AND ALL REPLICATIONS USING NO DISCOUNTING 
(DEVICE # 60) 
REGEST : OUTPUT FILE CONTAINING EACH PHASE ESTIMATE FOR EACH 
REPLICATION OF THE EXPONENTIAL REGRESSION ESTIMATE 
(DEVICE # 70) 


THE FOLLOWING IS A LIST OF KEY ARRAYS USED IN THE SIMULATION 


A ; MAIN WORKING ARRAY CONTAINS PROBABILITY OF SUCCESS FOR 
EACH FAILURE CAUSE, NUMBER OF TRIALS UNTIL FAILURE FOR 
EACH FAILURE CAUSE AND THE SYSTEM, CAUSE OF FAILURE, 
PHASE NUMBER, ADJUSTED NUMBER OF TRIALS AND ADJUSTED 
NUMBER OF FAILURES 
DIMENSION ( ((2%##CAUSES)+6) ,#FAILURES ) 

NFAPH : CONTAINS THE NUMBER OF FAILURES IN EACH PHASE 
DIMENSION (1,#PHASES) 

NFCAUS : BINARY ARRAY USED TO DETERMINE IF A FAILURE OCCURRED IN 
A PHASE 
DIMENSION ( 1,/FAILURE CAUSES) 

NTRIAL : CONTAINS THE NUMBER OF TRIALS SINCE LAST FAILURE OR 
DISCOUNTING FOR EACH FAILURE CAUSE 
DIMENSION ( 1,#FAILURE CAUSES ) 

TADJTP : CONTAINS THE TOTAL NUMBER OF ADJUSTED TRIALS IN A PHASE 
DIMENSION ( 1, # OF PHASES) 

PHREST : RECORDS THE PHASE ESTIMATE FOR EACH ESTIMATOR WITHIN A 
SINGLE REPLICATION 
DIMENSION (7,#PHASES) 
ROW 1 : WOODS WEIGHTED AVERAGE EST 
ROW 2 : MLE WITH DISCOUNTING 
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Saat = a - = - =. - = ae = = =@e ae a =a om =a a a a - T 4 - -. — - es se 
esededededisedededevededetoledededededestedededevedededesededlesedededektedeleeickiedekesdesededestesetideiedesetesesesese eee ee 


AREL 


CUMSF 


REG 


WREG1 


WREG2 
WREG3 


THE REMAINING ARRAYS ARE USED TO COMPUTE THE MEAN AND VARIANCE 
OF EACH ESTIMATE AT EACH PHASE. THEY ALL HAVE THE SAME DIMENSIONS 


ROW 3 = 
ROW 4&4 : 
ROW. 5: 
ROW on: 
ROW 7 : 


SINGLE PHASE MLE 


WOODS REGRESSION ESTIMATE (UNWEIGHTED) 
WEIGHTED REGRESSION ESTIMATE (METHOD 1) 
WEIGHTED REGRESSION ESTIMATE (METHOD 2) 
WEIGHTED REGRESSION ESTIMATE (METHOD 3) 
CONTAINS ACTUAL SYSTEM RELIABILITY IN EACH PHASE 


DIMENSION (1,#PHASES) 


CONTAINS THE NUMBER OF SUCCESS AND FAILURES FOR EACH 
FAILURE CAUSE (USED WITH WOODS WEIGHTED AVERAGE EST. ) 


DIMENSION (3,#FAILUE CAUSES) 


ROW 1 : 
ROW =2ee 
ROW 3 : 


NUMBER OF FAILURES 
NUMBER OF SUCCESSES 
ADJUSTED NUMBER OF SUCCESSES 


CONTAINS YJK VALUES UP TO 1000 
DIRE NSTON C15 1000) 


CONTAINS THE VARIANCE OF THE ESTIMATE AT EACH PHASE 


DIMENSION (1,##PHASES) 


CONTAINS THE RUNNING SUM OF 1/VAR(K) UP TO THE KTH PHASE 


DIMENSION (1,7#PHASES ) 


CONTAINS THE RUNNING SUM OF VAR(K) UP TO THE KTH PHASE 


DIMENSION (1,##PHASES) 


CONTAINS THE RUNNING SUM OF WEIGHTS UP TO THE KTH PHASE 


DIMENSION (1,/#PHASES) 


CONTAINS THE WEIGHT GIVEN TO EACH RELIABILITY ESTIMATE AT 


EACH PHASE. DIMENSION (3,#PHASES) 


ROW 1: WEIGHT = (1/VAR(K))/(SUM OF 1/VAR(K) UP TO 


ROW 2: WEIGHT = VAR(K)/(SUM OF VARCK) UP TO KTH PHASE) 


THE KTH PHASE) 


ROW 3: WEIGHT = USER INPUTTED VALUES. 


ARRAY USED TO 


ULES LOS (5 ,#PHASES) 


ROW 1 


ROW 2 : 
ROW 35% 
ROW 4 : 


ROW 5 


ARRAY USED TO COMPUTE THE WEIGHTED EXPONENTIAL REGRESSION 
BoUitATE USING METHOD 1 WEIGHTS. DIMENSION (4,/PHASES) 

ROW 1 : 
ROW 2 : 
ROWsS": 


ROW 4 


K BAR 
Y oBaAK 
Y BAR FOR THE PHASE 
B HAT 
A HAT 


X BAR(W) 
B HAT(W) 
A HAT(W) 
Y BAR(W) 


SAME AS ABOVE USING METHOD 2 WEIGHTS 
SAME AS ABOVE USING METHOD 3 WEIGHTS 


AND STRUCTURE 
DIMENSION (4,#PHASES ) 


MLEWD 
MLESP 


ROW 1 : 
ROW 2 : 
ROW 73. : 
ROW 4 : 


RUNNING SUM OF ESTIMATES 

RUNNING SUM OF SQUARED ESTIMATES 
MEAN OF THE ESTIMATES 

STANDARD DEVIATION OF THE ESTIMATES 


VALUES FOR THE MLE WITH DISCOUNTING 
VALUES FOR THE SINGLE PHASE MLE 


89 


COMPUTE THE EXPONENTIAL REGRESSION ESTIMATE 


* REGEST : VALUES FOR THE EXPONENTIAL REGRESSION ESTIMATE * 
* WREST1 : VALUES FOR THE WEIGHTED EXPONENTIAL REGRESSION ESTIMATE ¥ 


* USING METHOD 1 WEIGHTS. * 
* WREST2 : VALUES FOR THE WEIGHTED EXPONENTIAL REGRESSION ESTIMATE ¥* 
os USING METHOD 2 WEIGHTS. * 
* WREST3 : VALUES FOR THE WEIGHTED EXPONENTIAL REGRESSION ESTIMATE * 
a USING METHOD 3 WEIGHTS. * 


aetesesescoeacdedeseseuesdedesescakeseskocotcsesscokeskoescscseotseskscoktacscatacslkesescuese abacus sess aka aes seakseskslesesk okeseskesesescsteste 


C DEFINE AND DIMENSION VARIABLES 


PARAMETER (NR=50,NC=200) 

INTEGER REP ,DISOPT,FRELOP,LDI,ALD 

REAL*4 MIN 

REAL*8 DSEED,MLESP,MLEWD,EUL,SUM1,SUMS1,SUM2,SUMS2,SUM3 
C,SUMS3 

DIMENSION NFAPH(NR),A(NR,NC) ,NFCAUS(NR) ,NTRIAL(NR) , PHREST(7,NR), 
CMLEWD(4,NR) ,MLESP(4,NR) ,REGEST(4,NR) , AREL(NR) , YJK( 1000), 
CREG(5,NR) , TADJTP( 1000) , VAR(NR) , SVAR1(NR) ,W(3,NR), 
CWREG1(4,NR) ,WREST1(4,NR),WREST2(4,NR) ,WREST3(4,NR), 
CSVAR2(NR) , WREG2(4,NR) ,WREG3(4,NR) , SUMW(NR) ,EST(4,NR) , CUMSF( 3, NR) 


C READ IN THE NUMBER OF CAUSES TO BE USED ( NCAUSE ) AND THE NUMBER 
© OF PHASES GNPHASE IN tee Sd 


READ(10,*) NCAUSE 
READ(10,%*) NPHASE 


C CHECK IF FIXED RELIABILITY OPTION IS DESIRED. FIX EULER'S NUMBER. 


READ(10,*) FRELOP 
EUL = 0.5772156648 


CREATE VARIABLES FOR THE ROW INDICES OF THE WORKING MATRIX ( A ) 
IPHASE: PHASE 
ISYSPR: ACTUAL COMPONENT RELIABILITY 
INTR: NUMBER OF TRIALS UP TO AND INCLUDING FAILURE 
IFAILC: CAUSE OF THE FAILURE 
IADJF: ADJUSTED NUMBER OF FAILURES ED 
AFTER DISCOUNTING HAS BEEN APPLIED 
IADJT: ADJUSTED NUMBER OF TRIALS AFTER DISCOUNTING HAS BEEN APPLIED 
IYJK: YJK COMPUTED ON THE ADJUSTED NUMBER OF TRIALS 


SP ik SP SP ek SP SE a 


IPHASE 
ISiSPR 
INTR = 
IFAILC 


( 2*"NCAUSE)+1 
IPHASE +1 
SYSPR + 1 
INTR + 1 
IADJF = IFAILC + 1 
IADJT = IADJF + 1 
IYJK = IADJT + 1 


lew i 
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C READ IN THE NUMBER OF FAILURES IN EACH PHASE ( NFAPH(I) ) AND 
C COMPUTE THE TOTAL NUMBER OF FAILURES IN THE TEST ( NFAIL ) 


C 
C 


C 
C 


NFAIL = 0 
DO 10 I=1,NPHASE 
READ(10,*) NFAPH(CI) 
NFAIL = NFAIL + NFAPH(1) 
CONTINUE 


INPUT THE PROBABILITY OF SUCCESS IN A SINGLE TRIAL FOR EACH CAUSE 
IN EACH PHASE IF FRELOP EQUALS ONE. 


IF (FRELOP .EQ. 1) THEN 
DO 15 I=1,NCAUSE 
L =1 
DO 16 J=1,NPHASE 
READ(10,*) QQ 
DO 17 K = L,L+NFAPH(J)-1 


ACI, K) = QQ 
CONTINUE 
L = L+NFAPH(J) 
CONTINUE 
CONTINUE 
ELSE 


INPUT THE PROBABILITY OF SUCCESS IN A SINGLE TRIAL FOR EACH CAUSE 
PNSTHE FIRST PHASE If FRELOP EQUALS ZERO. 


DO 20 I=1,NCAUSE 
READ(10,*) A(I,1) 
CONTINUE 


ENDIF 


INPUT THE REMAINING VARIABLES , THE NUMBER OF SUCCESSFUL TRIALS 
BEPORE@A DISCOUNT IS APPLIED (N); THE DISCOUNT FACTOR (R); THE SEED 
FOR THE RANDOM NUMBER GENERATOR, GGUBFS, (DSEED); RELIABILITY 
GROWTH FRACTION (FRIMP); TRIGGER FOR PRINTING INTERMEDIATE OUTPUT 
(IOPT) 
ARIGGERS a SAVING EACH ESTIMATE AT EACH PHASE FOR EACH ESTIMATOR 

LOPy | WOODS WEIGHTED AVERAGE MODEL 

LORD Mor With DISCOUNTING 

ROETS eo INGLE PHASE MLE 

IOPT4 : EXPONENTIAL REGRESSION MODEL 

Ee veoh nD EXPONENITAL REGRESSION MODEL (METHOD 1) 

IOPT6 : WEIGHTED EXPONENTIAL REGRESSION MODEL (METHOD 2) 

TORT? WEIGHTED EXPONENTIAL REGRESSION MODEL (METHOD 3) 
DISCOUNTING OPTION TRIGGER (DISOPT); LLOYD FAILURE DISCOUNTING 
PARAMETER (GAMMA); LLOYD DISCOUNT INTERVAL 


READ(10,*) N 
READ(10,%) R 
READ(10,**) DSEED 
READ(10,*) FRIMP 
READ(10,**) NREP 
READ(10,**) IOPT 


a 


READ(10,*) IOPT1 
READ(10,*) IOPT2 
READ(10,**) IOPT3 
READ(10,*) IOPT4 
READ(10,*) IOPTS 
READ(10,*) IOPT6 
READ(10,*) IOPT7 
READ(10,*) DISOPT 
READ(10,*) GAMA 
READ(10,%) LD! 


G=-<<=5 READ IN THE USER INPUTTED WEIGHTS FOR COMPUTATION OF METHODS: 
DO 25 I=1,NPHASE 
READ(10,*) W(C3,1) 
25 CONTINUE 
ANREP = NREP 
DSEED1 = DSEED 


C INITIALIZE THE ARRAYS USED TO COMPUTE THE MEAN AND STANDARD DEVIATION 
C OF EACH ESTIMATOR 


DO 30 J=1,NPHASE 


DO 20ers 
IC py =, © 
MLEWD(I,J) = 0.0 
MLESP(I,J) = 0.0 
REGEST(I,J) = 0.0 
WREST1(I,J) = 0.0 
WREST2(1I,J) = 0.0 
WREST3(I,J) = 0.0 

30 CONTINUE 
DO 31 J=1,NPHASE 

DO 31 I=1,6 

PHREST(I,J) = 0.0 


31 CONTINUE 


C COMPUTE AND STORE THE YJK VALUESSUP 10° 10¢¢@ 


TRG) =. 0 
DO 40 1=1,999 
YIKCI41) = -YIKC Die or 
40 CONTINUE 


C COMPUTE AND STORE K BAR FOR THE EXPONENTIAL REGRESSION MODEL 


SUM = 0.0 
DO 50 I=1,NPHASE 
SUN = SUM: tal 
REGCiaje= SUNT 
a0 CONTINUE 


C MAJOR REPETITION OF THE SIMULATIONS EOOE 
DO 500 REP=1,NREP 


C INITIALIZE FAILURE CAUSE VECTOR (NFCAUS) AND (CUMSF) 
CeeCOtTrUTE THE INITIAL GrsTEN RELIABILITY 


REL = l. 
DO 60 I=1,NCAUSE 
NFCAUS(I) = 0 
REL = REL * A(T, 1) 
DO 60 J=1,3 
CUMSF(J,I) = 0 
60 CONTINUE 


C INITIALIZE COLUMN (FAILURE # ) COUNTER FOR THE WORKING ARRAY (A) 
J=1 


C LOOP TO COMPUTE THE NUMBER OF TRIALS UP TO AND INCLUDING FAILURE 
C AND THE CAUSE OF FAILURE FOR EACH FAILURE IN EACH PHASE 


DO 130 K=1,NPHASE 


C SKIP ACTUAL COMPCNENT RELIABILITY COMPUTATION AFTER FIRST REP 
C AND FOR FIRST FAILURE 


IF(J.EQ.1) GOTO 75 
IF(REP.GT.1) GOTO 75 
REL = 1. 


GC IF FIXED RELIABILITY OPTION IS SELECTED THEN PHASE RELIABILITIES 
C ARE COMPUTED AS FOLLOWS 


IF (FRELOP .EQ. 1) THEN 
DO 65 I=1,NCAUSE 
REL = REL*A(I,J) 
NFCAUS(I) = 0 
65 CONTINUE 
ELSE 


C COMPUTE NEW ACTUAL RELIABILITY FOR THE COMPONENT IN PHASE K 
DO 70 I=1,NCAUSE 


C INCREASE CAUSE PR(SUCCESS) IF IT CAUSED FAILURE IN THE PREVIOUS PHASE 
GC COMPUTE NEXT PHASE RELIABILITY AND REINITIALIZE NFCAUS (NOT USED IF 
fee DXED PHASE RELIABILITY OPTION IS SELECTED). 


TEGNFCAUSC1])-~ EG. 1) THEN 
AiO je—eaclaCo=1)) + (C1. = AC(I,(J-1)))*FRIMP) 
Pie vECNPCAUS(]).NE. 1) THEN 
Cl yl) = velco) 
ELSE 
ENDIF 
REL = REL*ACI,J) 
NFCAUS(I) = 0 
70 CONTINUE 


ENDIF 
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15 x de | 
TRTOT = 0.0 


C COMPUTE THE NUMBER OF TRIALS UP TO AND INCLUDING FAILURE AND THE 
C CAUSE OF FAILURE FOR EACH FAILURE IN THE PHASE 


DO 120 L=1,NFAPH(K) 
IF(REP.GT.1) GOTO 90 
IF(J1.EQ.1) GOTO 85 
IF (FRELOP .EQ. 1) GOTO 85 
DO 80 I=1,NCAUSE 

AC(I,J) = ACT,(J-1)) 


80 CONTINUE 
85 ACISYSPR,J) = REL 

AC(IPHASE,J) = K 
90 MIN = 7. 2E75 


DO 110 I=1,NCAUSE 
C ASSIGN # TRIALS FOR CAUSES WITH PR(SUCCESS) = 0 OR i 


IF(AC(1I,J).GE. 1.) THEN 
A((I+NCAUSE) ,J) = 7. 2E75 
GOTO 100 

ELSEIF(A(I,J).EQ.0.) THEN 
A((I+NCAUSE) ,J) = 1. 


GOTO 100 
ELSE 
ENDIF 
C CONVERT UNIFORM (0,1) RANDOM VARIABLE TO GEOMETRIC (# TRIALS UNTIL 
C FAILURE ) FOR EACH FAILURE CAUSE. RECORD THE MIN # TRIALS FOR THE 
C CAUSES AS THE SYSTEM # TRIALS UP TO AND INCLUDING FAILURE AND 
C RECORD THE FAILURE CAUSE 
ACC I+NCAUSE) ,J) = INT(1.+(LOG(GGUBFS(DSEED) )/LOG(A(I,J)))) 
100 IF(A(( I+NCAUSE) ,J). LE. MIN) THEN 
MIN = A((I+NCAUSE),J) 
IMIN = I 
ELSE 
ENDIF 
116 CONTINUE 


ACIFAILC,J) = IMIN 
NFCAUSCIMIN) = 1 


C COMPUTE THE TOTAL # OF TRIALS FOR THE MLE SINGLE PHASE ESTIMATE AND 
C INCREMENT FAILURE # COUNTERS 


AC INTR,J) = MIN 
TRTOT = TRTOT + ACINTR,J) 
SS ee 
Shi gpl 

120 CONTINUE 
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C COMPUTE THE MLE ESTIMATE OF COMPONENT RELIABILITY FOR THIS PHASE AND 
C COMPUTE THE RUNNING SUM OF ESTIMATES AND THE SUM OF ESTIMATES SQUARED 
CG FOR COMPUTATION OF THE MEAN AND STANDARD DEVIATION OF THE ESTIMATE 


PHREST(3,K) = (TRTOT - NFAPH(K))/TRTOT 
MLESP(1,K) = MLESP(1,K) + PHREST(3,K) 
MLESP(2,K) = MLESP(2,K) + (PHREST(3,K)7*2) 


130 CONTINUE 


CG INITIALIZE THE ADJUSTED NUMBER OF FAILURES TO 1 AND THE COUNT OF THE 
C NUMBER OF TRIALS SINCE FAILURE OR DISCOUNTING (NTRIALS(I) ) TO 0 
C IN PREPARATION FOR THE DISCOUNTING ROUTINE 


DO 140 J=1,NFAIL 
ACIADJF,J) = 1. 
140 CONTINUE 


DO 150 I=1,NCAUSE 
NTRIAL(I) = 0 
150 CONTINUE 


C INITIALIZE THE VARIANCE,VAR,THE SUM OF 1/VAR,SVAR,AND THE WEIGHTS,VW. 
DO 155 I=1,NPHASE 
SVAR1(I) 

SVAR2(I) 

SUMW(I) 

VAR(I) 

W(1,1) 

W(2,I) 


O27 72) =) 


155 CONTINUE 


DISCOUNTING ROUTINE REVIEWS ALL PAST FAILURES AND CAUSES TO DATE 
AND DETERMINES IF THE DISCOUNTING CONDITIONS HAVE BEEN MET. COMPUTES 
THE ADJUSTED FAILURES, THE ADJUSTED # OF TRIALS AND YJK 


(oy 1 a a ep) 


= 6 
DO 300 K=1,NPHASE 
DO 200 L=1,NFAPH(K) 
J=J+til 


C UPDATES THE NUMBER OF TRIALS SINCE FAILURE OR DISCOUNTING FOR EACH 
C FAILURE CAUSE 


ICAUSE = INT(A(IFAILC,J)+. 5) 
DO 160 I=1,NCAUSE 
IFC ICAUSE. EQ. I) THEN 
NTRIAL(I) = 0 
ELSEIF(ICAUSE. NE. I) THEN 
NTRIAL(I) = NTRIAL(I) + INT(ACINTR,J)+. 5) 
ELSE 
ENDIF 
160 CONTINUE 


200 CONTINUE 


C CHOOSE DISCOUNTING METHOD TO BE USED 
[PCUISOPTANES 2) GOTO Sikso 
C PERFORM LLOYD'S FAILURE DISCOUNTING METHOD 


DO 170 I=1,J 
Tl = INTCACIFAILC,1)+.5) 
IF(NTRIAL(I1).EQ.0) THEN 
ACIADJF,I) = 1.0 
elouiay agi 
ELSE 
ENDIF 


C THIS IS THE MODIFIED LLOYD METHOD USING A DISCOUNT INTERVAL. THE 
GC ORIGINAL DISCOUNT METHOD MAY BE EMPLOYED BY SETTING LDI TO ONE. 


ALD = INT(NTRIAL(1I1)/LDI) 
IF(ALD .EQ. 0) THEN 
ACIADJF,I) = 1.0 
GO ie 7s 
ELSE 
ACIADJF,I) = 1.0 - ((€1. -GAMA)**(1.0/ALD) ) 
ENDIF 


170 CONTINUE 
GOTO 210 


GC PERFORMS STRAIGHT PERCENT FAILURE DISCOUNTING AND 
C COMPUTES THE ADJUSTED # OF FAILURES 


180 DO 190 I=1,J 
Il = INT(ACIFAILC,I)+. 5) 
IF(NTRIAL(I1).EQ.0) THEN 
ACTADIr ye 
ELSEIF(NTRIAL(11).GE.N) THEN 
ACIADJF,I) = ACIADJF,1I)*((1. -R)**(NTRIAL(1I1)/N)) 
ELSE 
ENDIF 
190 CONTINUE 


C ADJUSTS THE # TRIALS SINCE FAILURE OR DISCOUNTING FOR THOSE CAUSES 
GC THAT HAVE MET OR SURPASSED THE DISCOUNTING THRESHOLD 
C FOR THE STRAIGHT PERCENT DISCOUNTING METHOD 


DO 205 I=1,NCAUSE 
IF(NTRIALCI).GE.N) NTRIALCI) = MODCNTRIALCT) ,N) 
205 CONTINUE 
210 TADJT = 0.0 
TYIK =50. 0 
TRY 0.0 
Kl = 0 


DOZ15, 12=1,53 


DO 215 I=1,NCAUSE 
CUMSF(12,I) = 0 
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PAS CONTINUE 


C COMPUTES THE ADJUSTED # OF TRIALS FROM THE ADJUSTED # OF FAILURES 
C AND COMPUTES THE SUM OF THE ADJUSTED # OF TRIALS FOR ESTIMATE COMP. 


PREL = 0.0 
LTRIAL = 0 


C IF ADJUSTED FAILURES ARE APPROACHING 0 THEN ADJUSTED TRIALS MUST 
teecr PRE-SET. 


DO 240 I=1,J 


IF(ACIADJF,I) .LE. .0000001) THEN 
AC IADJF,I) = .0000001 
ENDIF 


AC IADJT,I) = ACINTR,I)/ACIADJF,1) 
ion = TADIn +eACTADIT, I> 


C COMPUTE YJK FROM THE ADJUSTED # OF TRIALS AND STORE THE SUM FOR 
C ESTIMATE COMPUTATION, USE ARRAY FOR # TRIALS < 1000 AND APPROX. FOR 
C VALUES > 1000 


Nl = NINT(AC(IADJT,1I)) 
Pe Gio LOOO a HEN 
AOLYIK, 1) = YIKCNI) 
ESEIFCNI. GT. 1000) THEN 
X=N1 
Q=12*xK 
T=X+1 
S=A+2 


eves 1) -CRUESCEOG( %))=C1/ (252 Gl (G*T) = (1/(0"T*S))) 


ELSE 
SNE 


C DETERMINE IF A PHASE BOUNDARY HAS BEEN REACHED TO BEGIN ESTIMATE 
C COMPUTATION 


Leite. 1) GOTO.225 
IF( AC IPHASE ,I).NE. AC IPHASE,(I-1))) THEN 


C COMPUTE THE WOODS WEIGHTED AVERAGE ESTIMATE 


MAX = 0 
Kl = Kl + 1 


C DETERMINE THE FAILURE CAUSE WITH THE LARGEST # OF FAILURES 


DO 220 I1=1,NCAUSE 
IF(CUMSF(1,I1).GT.MAX) THEN 
MAX = CUMSF(1,11) 

ICOL = I1 
ELSE 


oF 


END TE 
220 CONTINUE 


C COMPUTE YJK VALUE FOR THE CURRENT PHASE ESTIMATE 


IF(CUMSF(1,ICOL). LE. 1000) THEN 
AHATL = YJKCCUMSF(1,ICOL)) 
ELSEIF(CUMSF(1,ICOL).GT. 1000) THEN 
X = CUMSF(1,ICOL) 
Q= te ee 
=X+1 
S=Xt2 


AHATL=(EUL+( LOG(X) )+01/(2*%) =(1/(Q"1) Gi Goes 


ELSE 
ENDIF 
IX = CUMSF(1,1COL) = CUMSE@3 Seay 
IF(IX. LE. 1000) THEN 

AHATU = YJK(C IX) 
ELSEIFC IX. GT. 1000) THEN 

X = IX 

Q=1 pa 

T=X+1 

S=Xt2 


AHATU=( EULt( LOG(X) )+(1/( 2" 2) ist TG ey ead er) 


ELSE 
ENDIF 


C COMPUTE CURRENT PHASE RELIABILITY ESTIMATE 
AHAT = AHATU - AHATL 


CREL 120) =VERPG@-4sHsT) 
®, = CUMSFCI], ICOL) > CURSEC3, Icom) 


C COMPUTE AND STORE THE WOODS WEIGHTED AVERAGE ESTIMATE 


PREL = ((LTRIAL*PREL)/X) + (((X-LTRIAL)*CREL) /X) 
LTRIAL = CUMSF(1,ICOL) + CUMSF(3,ICOL) 


COMPUTE THE PHASE AND GLOBAL AVERAGE FOR YJK USED IN THE EXPONENTIAL 
REGRESSION ESTIMATES ARE 


oa 


REG(2,K1) 
REG(3,K1) 
TPYJK = 0. 


= TYJK/(I-1) 

= TPYJK/NFAPHC(K1) 
0 

ENDIE 


C COMPUTE THE NUMBER OF FAILURES AND SUCCESSES FOR EACH FAILURE CAUSE 
C USED IN THE WOODS WEIGHTED AVERAGE ESTIMATE 


225 ICAUSE = INT(A(IFAILC,I)+. 5) 
DO 230 I1=1,NCAUSE 
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CUMSF(2,11) = CUMSF(2,11) + INT(ACINTR,I) + .5) 
CUMSF(3,I1} = CUMSF(3,I1) + Nl 
230 CONTINUE 
CUMSF(1,ICAUSE) = CUMSF(1,ICAUSE) + 1 
CUMSF(2,ICAUSE) = CUMSF(2,ICAUSE) - 1 
CUMSF(3,ICAUSE) = CUMSF(3,ICAUSE) - 1 
TPYJK = TPYJK + ACIYJK,1) 
TYJK = TYJK + A(IYJK,I) 
240 CONTINUE 


C REPEAT COMPUTATIONS FOR THE WOODS WEIGHTED AVERAGE ESTIMATE FOR THE 
C FINAL PHASE 


MAX = 0 

K1 = K1 + 1 

DO 245 I1=1,NCAUSE 
IF(CUMSF(1,11).GT.MAX) THEN 
MAX = CUMSF(1,I1) 
ICOL = 11 
ELSE 
ENDIF 

245 CONTINUE 


IPCCUMSEC 1, ICOiny LE. 1000) THEN 
AHATL = YJK(CUMSF(1,ICOL)) 
ELSEIFCCUMSF(1,ICOL).GT. 1000) THEN 
X = CUMSF(1,ICOL) 
Q=12*xX 
T=X+1 
S=K+2 


eae ee EUR CROC) (2*X) )-C1/CQ*T) )-(1/CQ*T*S) ) } 


alee 
ENDIF 
Pe eecuUronC leon) + CUMSFC3,ICOL) 
PGi be. 1000) THEN 

AHATU = YJKCIX) 
ELSEIFC IX. GT. 1000) THEN 

X = IX 

Q=12*X 

=X+1 

S=X+2 


AHATU=( EUL+( LOG(X) )+01/C2*X) ) -C1/CQ*T) ) -(1/CQ*T*S ))) 


ELSE 
ENDIF 

AHAT = AHATU - AHATL 
Coie or= PPC AHA) 


ReaeeUlon( 1, 1COL) + CUMSFC3,1COL) 


PREL = ((LTRIAL*PREL)/X) + (((X-LTRIAL)*CREL) /X) 
LGRieie= CUMSEC1,1COL) + CUNSF(3,1COL) 
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REG(2,K1) = TYJk ep 
REG(S,K1) = TPYJK/NFAPHCK® 


PHREST(1,K) = PREL 
C COMPUTE THE MLE ESTIMATE OF PHASE RELIABILITY USING DISCOUNTING 
PHREST(2,K) = (TADJT - J)/TADJT 
C COMPUTE THE EXPONENTIAL REGRESSION ESTIMATE BEGINNING WITH B HAT 
SUM = 0.0 
SUMS = 0.0 
IF (K. BO. 1) GOnem252 
be) OS) ae 
SUM = SUM + ((I-REG(1,K))*REG(3,1)) 
SUMS = SUMS + ((I-REG(1,K))**2) 
250 CONTINUE 
REG(4,K) = SUM/SUMS 
C COMPUTE A HAT 


REG(5 ,K) 


REG(2,K) - (REG(4,K)*REG(1,K)) 
C COMPUTE AND STORE THE EXPONENTIAL REGRESSION ESTIMATE 


PHREST(4,K) = 1.0 - EXP(-(REG(5,K) + (REG(4,K)*K))) 
IF( PHREST(4,K). LT. 0.0) PHREST(4,K)=0. 0 
GOTO 255 

252 PHREST(4,K) = 1.0 - EXP(-REG(3,1)) 
IF( PHREST(4,K). LT. 0.0) PHREST(4,K)=0. 0 


STORE THE RUNNING SUM OF THE ESTIMATES FOR THE CURRENT PHASE AND THE 
RUNNING SUM OF THE ESTIMATES SQUARED FOR COMPUTATION OF THE MEAN AND 
STANDARD DEVIATION OF EACH ESTIMATE FOR EACH RELIABILITY GROWTH 
MODEL 


Ca Cae) Ca 


255 EST(1,K) = EST(1,K) + PHREST(1,K) 
EST(2,.K) = EST(2.K) + (PHRRenC ero 
MLEWD(1,K) = MLEWD(1,K) + PHREST(2,K) 
MLEWD(2,K) = MLEWD(2,K) + (PHREST(2,K)**2) 
REGEST(1,K) = REGEST(1,K) + PHREST(4,K) 
REGEST(2,K) = REGEST(2,K) + (PHREST(4,K)**2) 


C COMPUTE THE VARIANCES OF THE UNWEIGHTED EXPONENTIAL REGRESSION 
C ESTIMATES AND STORE THE RUNNING SUM OF THE VARIANCES FOR USE IN THE 
C WEIGHTED REGRESSION MODEL. 


VAR(K) = (NFAPH(K)*PHREST(1,K))/((1. -PHREST(1,K))*%2) 
IF ( VAR(K) .LT. .0000001)THEN 

VAR(K) = .0000001 

END IF 

IP UGK SEOs alpeGO Ton256 

SVARI(K) = SVARI(K-1) + (1. /VAR(K)) 
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250 


SVAR2(K) = SVAR2(K-1) + VAR(K) 
SUMW(K) = SUMW(K-1) + W(3,K) 
CO Ge259 


SVAR1(K) 
SVAR2(K) 


SUMW(K) 


C COMPUTE THE 
259 W(1,K) = 


W(2,K) = 


C COMPUTE AND 


261 


SUMX1 
SUP a 
SUMX2 
SUMi2 
SUMX3 
SUMY3 


DOwZCrs I 


SUMX1 
SUMY1 
SUMX2 
SUMY2 
SUMX3 
Shy 3 


CONTINUE 


(1. /VARC1)) 


VAR( 1) 
W(3,1) 


WEIGHTS FOR EACH PHASE 

(1. /VAR(K))/SVAR1(K) 

VAR(K) /SVAR2(K) 

STORE X-BAR(W) AND Y-BAR(W) 


OO eo oc 2] @& 


WREG1(1,K) 
WREG1(4,K) 
WREG2(1,K) 
WREG2(4,K) 
WREG3(1,K) 
WREG3(4,K) 


Ts 

SUMX1 
SUMY1 
SUMX2 
SUMY2 
SUMX3 
SUNY3 


+ (1. /VAR(K))*I 

+ (1. /VAR(K))*REG(3,1) 
+ VAR(K)*I 
VAR(K)*REG( 3,1) 
W(3,K)*I 


+ 
+ 
+ W(3,K)*REG(3,I) 


SUMX1/SVAR1(K) 
SUMY1/SVAR1(K) 
SUMX2/SVAR2(K) 
SUMY2/SVAR2(K) 
SUMX3 /SUMW(K) 
SUMY3/SUMW(K) 


C COMPUTE AND STORE B-HAT(W) 


27) 


SUM1 
SUMS1 
SUM2 
SUMS2 
SUM3 
SUMS3 


Chee. 


BOs Il 


SUM1 


SUMS1 


SUM2 


SUMS2 


SUM3 


SUMS3 
CONTINUE 


0. 
O 
QO. 
OF 
OQ. 
Q. 
OF 


WREG1(2 ,K) 
WREG2(2,K) 


SSS So o-s 


LP oO nO. 27 2 


Aves 


SUM1 + (W(1,K)*( I-WREG1(1,K))*REG(3,I)) 
SUMS1 + ( W(1,K)*((I-WREG1(1,K))***2)) 
SUM2 + (W(2,K)**(I-WREG2(1,K))*REG(3,I)) 
SUMS2 + ( W(2,K)*(( I-WREG2(1,K))#**2)) 
SUM3 + (W(3,K)*( I-WREG3(1,K))*REG(3,1)) 
SUMS3 + ( W(3,K)*((I1-WREG3(1,K))#*2)) 


SUM1/SUMS1 
SUM2 /SUMS2 
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WREG3(2,K) = SUM3/SUMS3 


C COMPUTE AND STORE A-HAT(W) 


WREGI1(3,K) = WREGI(4,K) - CWREGI(2,K)“WREGiGe a 
WREG2(3,K) = WREG2(4,K) = CWREG2(2,K)*WREGZ a 
WREG3(3,K) = WREG3(4,K) - (WREG3(2,K)*WREG3(1,K)) 


C COMPUTE AND STORE THE WEIGHTED EXPONENTIAL REGRESSION ESTIMATE 


PORES Cae ks) 
PHRES PC6a 0) 
PHREST(7,K) 


1.0 - EXP(-(WREG1(3,K)+(WREG1(2,K)*K))) 
1.0 - EXP(-(WREG2(3,K)+(WREG2(2,K)*K))) 
1.0 - EXP(-(WREG3(3,K)+(WREG3(2,K)*K))) 


IF (PHRESTCS5,K) . LT.* 0.0) @PaRe sia) — oe 
IF (PHREST(6,K) . UT. 0: 0) SPHRE SiGe = Seven 
IF (PHREST(7,K) . LT. 0.0) FPHREST C7) =a 
GOe TO e275 

272 PHRESI(5,K) = 1.0 - EXPC=WREE ICG aD) 
IF (PHREST(5,K) .&T..-0.0) PHREST Cs os 
PHREST(6,K) = 1, 0 ="EXPC-WREG2 Ge 
IF (PHREST(C6;K) 2 UT. 0: 0) PHREST GG, aed 
PHREST(7,K) = 1.0 = EXP( -WheGsG. amp 
IF (PHREST(7,K) .LT. 0.0) PHRESI@ tie ee ae 


GC STORE THE RUNNING SUM OF THE WEIGHTED ESTIMATES FORK FHESCURRENT 

C PHASE AND THE RUNNING SUM OF THE ESTIMATES SQUARED FOR COMPUTATION OF 
C THE MEAN AND STANDARD DEVIATION OF EACH WEIGHTED ESTIMATE FOR THE 

C WEIGHTED EXPONENTIAL REGRESSION MODEL. 


975 WREST1(1,%) 
WREST1(2,K) 


= WREST1(1,K) + PHREST(5,K) 

= WREST1(2,K) + (PHREST(5,K)**2) 
WREST2(1,K) 
WREST2(2,K) 


WREST2(1,K) + PHREST(6,K) 
WREST2(2,K) + (PHREST(6,K)**2) 


WREST3(1,K) 
WREST3(2,K) 


WREST3(1,4) =SPHRES Tee) 
WREST3( 2,4) + CRHRES Wye) 2) 


C STORE THE ACTUAL PHASE RELIABILITY 
AREL(K) = AC ISYSPR,J) 


C PRINT INTERMEDIATE OUTPUT IF REQUESTED AND THE NUMBER OF REPETITIONS 
C185 NOY GREATER InAN@S 


TPCIORT: NE. 1) GOles 390 
EECKEP.-GT.5)..GOT0" 200 


WRITE( 20,1000) REP,K 
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1000 FORMAT(T16,'REPETITION NUMBER: ',I4,' PHASE NUMBER: ',I4) 
WRITE( 20,1010) AC ISYSPR,J) 

1010 FORMAT( 22X,'ACTUAL COMPONENT RELIABILITY: '‘,F7.5) 
WRITE( 20,1020) PHREST(1,K) 

1020 FORMAT(20X, ‘PREDICTED COMPONENT RELIABILITY: ',F7.5) 
WRITE( 20,1022) PHREST(2,K) 

1022 FORMAT(20X,'MLE ESTIMATE USING DISCOUNTING: ',F7.5) 
WRITE( 20,1025) PHREST(3,K) 

1025 FORMAT(18X,'MLE ESTIMATE OF PHASE RELIABILITY: ',F7.5) 
WRITE( 20,1027) PHREST(4,K) 

1027 FORMAT(14X, 'REGRESSION ESTIMATE OF PHASE RELIABILITY: ',F7.5) 
WRITE( 20,1028) PHREST(4,K) 


1028 FORMAT(14X,'WEIGHTED REG. ESTIMATE ( METHOD 1) =: ',F7.5) 
WRITE( 20,1026) PHREST(5,K) 

1026 FORMAT(14X,'WEIGHTED REG. ESTIMATE ( METHOD 2) : ',F7.5) 
WRITE(20,1029) PHREST(6,K) 

1029 FORMAT(14X,'WEIGHTED REG. ESTIMATE ( METHOD 3) : ',F7.5) 


WRITE( 20,1030) 
1030 FORMAT(' ',' 
DO 260 I=1,NCAUSE 
WRITE(20,1035)1I,A(1,J),AC(( I+NCAUSE) ,J) 
1035  FORMAT(12X,'CAUSE: ',13,' BRO SUCCESS): = F 7.6," ## TRIALS: ', 
Crlgee 
260 CONTINUE 
WRITE( 20,1036) 
1036 FORMAT(' ',' ) 
WRITE( 20,1040) 
ieaOuPOR enCax. PAIL # 93X, FAIL CAUSE’ ,3X,'# TRIALS’ ,3X,'ADJ # FAIL' ,3 
Gopobs + TRIALS: .7xX, YIK') 
Ome ae I—1.J 
Vee Oe 1050) 1 A [PADLG, 1),ACINDR,1),ACIADJF,1),ACIADJT,1I),ACIYJK 
Cal) 
o> 0 PORMATC 2X, 13,6X,F3.0,7X,F8.0,4X,F8.6,4X,F12. 0,3X,F11. 4) 
270 CONTINUE 
WRITE( 20, 1060) 
1060 FORMAT(' ',///) 


300 CONTINUE 
Peek Nisbeon OF THE 3 ESTIMATES TO THEIR APPROPRIATE OUTPUT FILE 
ter tt KEQUESTED 


IF(IOPT1.NE.1) GOTO 401 
400 WRITE(40,2000) (PHREST(1,I), I=1,NPHASE) 
401 IFCIOPT2.NE.1) GOTO 402 

WRITE(50,2000) (PHREST(2,1), I=1,NPHASE) 
402 IF(IOPT3.NE.1} GOTO 403 

WRITE(60,2000) (PHREST(3,1), I=1,NPHASE) 
403  IFCIOPT4.NE.1) GOTO 404 

WRITE( 70,2000) (PHREST(4,1I); I=1,NPHASE) 
404 IFCIOPT5.NE.1) GOTO 405 

WRITE( 16,2000) (PHREST(5,I), I=1,NPHASE) 
405 IF(IOPTS.NE.1) GOTO 406 

WRITE(39,2000) (PHREST(6,1), I=1,NPHASE) 
406 IFCIOPTS.NE.1) GOTO 500 

WRITE( 49,2000) (PHREST(7,1), I=1,NPHASE) 
2000 FORMAT(' ' ,30(F7. 6: 1X)) 


103 


500 CONTINUE 


C WORKING A MATRIX 
DO 4050 J = 1,(2*NCAUSE)+7 
WRITE(83,*) (ACJ,1), 1 = SINE ey 
4050 CONTINUE 


DO 4052 J = 1,NCAUSE 
WRITE(88,*) CREG(S,1)[>) = Fae) 
4052 CONTINUE 


GC UPON COMPLETION OF ALL REPETITIONS, COMPUTE THE MEAN AND STANDARD 
C DEVIATION OF EACH ESTIMATE FOR EACH PHASE SKIPPING COMPUTATIONS IF 
GC ONLY ONE REPETITION IS REQUIRED 


LE? (NREP. LE) COrOeeo 


DO 600 I=1,NPHASE 
EST(3,1) = EST(1,1)/XNREP 
MLEWD(3,1) = MLEWD(1,1)/XNREP 
MLESP(3,1) = MLESP(1,1)/XNREP 
REGEST(3,1) = REGEST(1,1)/XNREP 
WREST1(3,1) = WREST1(1,1)/XNREP 
WREST2(3,1) = WREST2(1,1)/XNREP 
WREST3(3,1) = WREST3(1,1)/XNREP 
EST(4,1) = SQRT((EST(2,1)-(XNREP*(EST(3,1)7**2)))/(XNREP-1)) 
MLEWD(4,1) = SQRT((MLEWD(2,1)-(XNREP**(MLEWD(3,1)2**2)))/CXNREP-1)) 
MLESP(4,1) = SQRT((MLESP( 2 ,1I)-(XNREP**(MLESP(3,1)*"2)))/CXNREP-1)) 
REGEST(4 , I )=SQRT( (REGEST(2 ,1)-(XNREP**( REGEST( 3, 1)**2)))/(XNREP-1)) 
) 
) 
) 


1D) 
WREST1(4, I)=SQRT( (WREST1(2,1)-(XNREP**(WREST1(3,1)**2)))/(XNREP-1)) 
WREST2(4, 1 )}=SQRT( (WREST2(2,1)-(XNREP**(WREST2(3,1)***2)))/(XNREP-1)) 
WREST3(4, I)=SQRT( (WREST3(2,1)-(XNREP**(WREST3(3,1)**2)))/(XNREP-1)) 
600 CONTINUE 


G PRINT THE FINAL OUTPUT TABLE 7O AVFILE 


601 WRITE(30,3000) 
3000 FORMAT('0' ,T47,'DISCRETE RELIABILITY GROWTH SIMULATION’ ) 
WRITE( 30, 3010) 
3010 FORMAT('-' ,T54,'MODEL PARAMETER SUMMARY’ ) 
WRITE( 30,3020) NCAUSE 
3020 FORMAT('0' ,T47,'NUMBER OF POSSIBLE FAILURE CAUSES ',I4) 
IF (FRELOP .EQ. 1) GOTO 4000 
WRITE( 30, 3030) 
3030 FORMAT('0' ,T38,'CAUSE NUMBER' ,T64,'SINGLE TRIAL PR( SUCCESS ) FOR 
CPHASE 1') 
DO 3050 M=1,NCAUSE 
WRITE( 30,3040) M,A(M,1) 
3040 FORMAT(' ‘,T43,12,T79,F8. 6) 
3050 CONTINUE 
WRITE( 30,3060) FRIMP 
3060 FORMAT('0' ,T37,'FRACTION CAUSE RELIABILITY IMPROVES AFTER FAILURE 
Cy, Feeo) 
5000 WRITE(30,3080* NPHASE 
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3080 FORMAT('-',T48,'NUMBER OF PHASES IN THE SIMULATION ' ,1I2) 

WRITE( 30, 3090) 
3090 FORMAT('0' ,T42,'PHASE NUMBER’ ,T59,'NUMBER OF FAILURES IN THE FIRST 

C PHASE') 

DO 3110 M=1,NPHASE 

WRITE( 30,3100) M,NFAPH(M) 
BOO eRORAIM sees .12,773,12) 
3110 CONTINUE 

WRITE( 30,3120) NFAIL 
3120 FORMAT('0',T51,'TOTAL NUMBER OF FAILURES ' ,I4) 

IF(DISOPT. EQ. 2) GO TO 3160 


WRITE( 30,3130) 
3130 FORMAT('-',T38,'DISCOUNTING PERFORMED USING THE CONSTANT FRACTION 
CMETHOD' ) 


WRITE( 30,3140) R 
3140 FORMAT('0',T44,'FRACTION EACH FAILURE IS DISCOUNTED ',F8. 6) 
WRITE( 30,3150) N 
3150 FORMAT(' ',T33,'NUMBER OF TRIALS AFTER A FAILURE BEFORE A DISCOUNT 
C IS APPLIED ',I4) 
GO TO 3190 
3160 WRITE( 30,3170) 
3170 FORMAT('-' ,T44,'DISCOUNTING PERFORMED USING THE LLOYD METHOD' ) 
WRITE( 30,3180) GAMA 
3180 FORMAT('0' ,T39,'PERCENT C.I. ( USED AS DISCOUNT FRACTION ) ',F8.6 
C) 
WRITE( 30,3185) LDI 
3185 FORMAT('0' ,T50,' LLOYD DISCOUNT INTERVAL: ‘ ,1I3) 
3190 WRITE( 30,3200) DSEED1 
3200 FORMAT('-',T46,' RANDOM NUMBER SEED USED ° ,F15. 2) 
WRITE( 30,3210) NREP 
3210 FORMAT('0' ,T37,' NUMBER OF REPETITIONS OF THE SIMULATION PERFORMED 
Orly) 
WRITE( 30,3220) 
3220 FORMAT('1' ,T61, 'ESTIMATOR: ') 
WRITE( 30,3230) 
3230 FORMAT('0' ,T48,'SINGLE PHASE MLE WITHOUT DISCOUNTING’ ) 
WRITE( 30,3240) 
3240 FORMAT('-' ,T60, 'MEAN' ,T83, ‘ESTIMATE’ ,T109,'95 %') 
WRITE( 30,3250) 
3250 FORMAT(’ ',T12,'PHASE NUMBER' ,T29,' ACTUAL RELIABILITY’ ,T52,' PREDIC 
CTED RELIABILITY' ,T78,' STANDARD DEVIATION’ ,T101, CONFIDENCE INTERVA 
Sl} 


Ee COMPUTE €.1. FOR SINGLE PHASE MLE 


DO 3270 M=1,NPHASE 
CI = (1. 96%*MLESP(4,M))/SQRT( XNREP) 
CIU = MLESP(3,M) + CI 
CIL = MLESP(3,M) - CI 
WRITE( 30,3260) M,AREL(M) ,MLESP(3,M),MLESP(4,M),CIL,CIU 
OOmm@RMaiK 0 ,017,12,934,F8.6,158,F8.6,182,F9.6,799,'( ',F8.6,' , ',F 
Clare) ) 
3270 CONTINUE 
WRITE( 30,3220) 
WRITE( 30,3280) 
3280 FORMAT('0',T42,'MAX LIKELIHOOD ESTIMATE USING DISCOUNTED FAILURES' 


C) 
WRIdn@e0 5240) 
WRETEC O0e3a 250) 


C COMPUTE C.I. FOR MLE WITH DISCOUNTING 


DO 3290 M=1,NPHASE 
CI = (1. 96*MLEWD(4 ,M) ) /SQRT( XNREP) 
CIU = MLEWD(3,M) + CI 
CIL = MLEWD(3,M) - CI 
WRITE( 30,3260) M,AREL(M) ,MLEWD(3,M) ,MLEWD(4,M) ,CIL,CIU 
3290 CONTINUE 
WRITE( 30,3220) 
WRITE( 30, 3300) 
3300 FORMAT('0' ,T38,' WEIGHTED AVERAGE ESTIMATE USING FAILURE DISCOUNTIN 
CG') 
WRITE( 30,3240) 
WRITE( 30, 3250) 


C COMPUTE C.I. FOR WOODS WEIGHTED AVERAGE ESTIMATES 


DO 3310 M=1,NPHASE 

CI = (1. 96*EST(4,M))/SQRT(XNREP) 

CIU = EST(3,M) + CI 

CIL = EST(3,M) - CI 

WRITE( 30,3260) M,AREL(M),EST(3,M),EST(4,M) ,CIL,CIU 
3310 CONTINUE 


WRITE( 30,3220) 
WRITE( 30 ,3320) 
3320 FORMAT('0' ,T43,'REGRESSION ESTIMATE USING DISCOUNTED FAILURES’ ) 
WRITE( 30,3240) 
WRITE( 30, 3250) 


C COMPUTE C.I. FOR EXPONENTIAL REGRESSION ESTIMATES 


DO-3350; M=1,NPHASE 

CI = (1. 96*REGEST(4,M))/SQRT( XNREP) 

C1U, ={REGEST C3, Mes 61 

Cit = REGESTCS.) =- cl 

WRITE(30 ,3260) M,AREL(M) ,REGEST(3 ,M) ,REGESI( 4,0) (Glin ao 
5330 CONTEUE 


C----- WEIGHTED REGRESSION (METHOD 1) 
WRITE( 30, 3220) 
WRITE( 30,3321) 
3321 FORMAT('0' ,T43,'WEIGHTED REGRESSSION ESTIMATE (METHOD 1 )') 


WRITE( 30, 3240) 
WRITE( 30,3250) 


C COMPUTE C.J]. FOR WEIGHTED EXPONENTIAL REGRESSION ESTIMATES 


DO 3331 M=1,NPHASE 
CI = (1. 96*WREST1(4,M))/SQRT(XNREP) 
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CIU = WREST1(3,M) + CI 

CIL = WREST1(3,M) - CI 

WRITE( 30,3260) M,AREL(M) ,WREST1(3,M) ,WREST1(4,M) ,CIL,CIU 
3331 CONTINUE 


C----- WEIGHTED REGRESSION (METHOD 2) 
WRITE( 30,3220) 
WRITE( 30, 3322) 
3322 FORMAT('0',T43,'WEIGHTED REGRESSION ESTIMATE (METHOD 2) ') 


WRITE( 30,3240) 
WRITE(30, 3250) 


C COMPUTE C.I. FOR WEIGHTED EXPONENTIAL REGRESSION ESTIMATES 


DO 3332 M=1,NPHASE 

CI = (1. 96*WREST2(4,M))/SQRT(XNREP) 

CIU = WREST2(3,M) + CI 

CIL = WREST2(3,M) - CI 

WRITE( 30,3260) M,AREL(M) ,WREST2(3,M) ,WREST2(4,M) ,CIL,CIU 
3332 CONTINUE 


C----- WEIGHTED REGRESSION (METHOD 3) 
WRiMEGsOs220) 
WRITE( 30, 3323) 
3323 FORMAT('0',T43, ‘WEIGHTED REGRESSION ESTIMATE (METHOD 3) ') 


WRITE( 30,3240) 
WRITE( 30,3250) 


C COMPUTE C.I. FOR WEIGHTED EXPONENTIAL REGRESSION ESTIMATES 


DO 3333 M=1,NPHASE 

Ciel oo: WREST 3(4,M)) /SQRTCXNREP) 

Ciie— Whee l3(3,M) + Cl 

Cit = Wreel3(3,M) = cl 

Wh CSO S260) M SAREL(M) ,WREST3(3),M)QWRESI3(4,M),CIL,CIU 
3333 CONTINUE 


WRITE( 30, 3340) 

3340 FORMAT( ‘1° ,T59, 'RECAPITULATION' //) 
WRITE( 30, 3350) 

3350 FORMAT('-' ,T3, PHASE’ ,T11, 'ACTUAL' ,T28, 'MEAN' ,T38, 'EST' ,T50, 'MEAN' 
mic ol 77s MEAN (182, EST .T92, MEAN ,T102,'EST' ,T112 
Ce MbANs, 1122, EST ) 
WRITE( 30,3360) 

SeCcO HORNAT(C ‘Till, RELIAB’ ,128,'WGT' ,T38,'STD' ,T50,'MLE' ,T62,'STD’ ,T7 
Tait wholy sib ,092, REG ,1T102, STD’ 57112,'WT' ,T122,'STD') 
WRITE( 30, 3370) 

e570) FORMAT( ',1T28, AVG ,1T38,'DEV' ,T50,'W/D' ,T62,'DEV' ,T72 
Pete woe DEV 6192, EST’ ,T102, DEV’ ,T112,'REG' ,T122, 'DEV') 
WRITE( 30,3375) 

Se75 FORMATO ,128, EST ,T112,'(1)'/) 
DO 650 J=1,NPHASE 

Nie sor 33c0 lL AREIG1),bo0(G.1),ES1(4,1),MLEWD(3,1) ,MLEWD(4,1), 

CMLESP(3,1),MLESP(4,1),REGEST(3,1),REGEST(4,1),WREST1(3,I1), 
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CWREST1(4,1) 
3380 FORMAT('0' ,T4,13,T11,F7.6,T26,F7. 6,136,F/. 6,150,F7.6,102,57 soma 
CF7.6,T82,F7.6,192,F7.6,T102,F7. 6,1112,F7. 6, oe 


650 CONTINUE 
DSSS PAGE 2 OF RECAPITULATION: 


WRITE( 30,3381) 
3381 FORMAT('1',T59,'RECAPITULATION CONT. ‘'//) 


WRITE( 30, 3382) 

3382 FORMAT('-',T3,'PHASE' 
CiNEAN alos. ESD — hyo. 
WRITE( 30, 3383) 

3383 FORMATC(' ',T1i, ‘RELIAB' 
C SUD 077. MEDHOD Tees 
WRITE( 30, 3384) 

3384 FORMAT(' ',T28,'EST' ,T35, DEVIATION’ ,1T53, ONE ,7T60, DEVIATIGRs 
C,T78,'TWO' ,T85, DEVIATION’ ,T103, ‘THREE , 1110) DEVieAnnon 
DO 651 I=1,NPHASE 
WRITE( 30,3385) I,AREL(1I),REGEST(3,1),REGEST(4,1) ,WREST1(3,I1) 
C,WREST1(4,1) ,WREST2(3,1),WREST2(4,1) ,WREST3( 3,1) ,WREST3(4,1) 

3385 FORMAT( 0° ,T4,13,T11,F7.6,1T26,F7.6,136,F/.6,151,F7.06, loler ae 
C,176.F7.6,186.F 726.1101" Py 6-0) lee 


T11,'ACTUAL' ,T28, MEAN’ (1385 Eon 


MEAN‘ ,T88, EST’ ,T103, ‘MEAN’ , T1134)" Boum 


; 
T28,'REG' ,T38,'STD' ,T53, ‘METHOD |toom 


‘stp! 1103, 'METHOD' ,T113,'STD") 


WRITE (53,3400) MLEWD(4,1) 
WRITE. (54,3400) MLEWD( 3.4) 


WRITE (63,3400) MLESP(4,1) 
WRITE (64,3400) MLESP(3,1) 


WRITE (73,3400) REGEST(4,I) 
WRITE (74,3400) REGEST(3,1) 


WRITE(17,3400) WREST1(4,1) 
WRITE(18,3400) WREST1(3,1) 
WRITE(93,3400) WREST2(4,1) 
WRITE( 94,3400) WREST2(3,I) 
WRITE(95,3400) WREST3(4,1) 
WRITE(96,3400) WREST3(3,1) 


WRITE (84,3400) AREL(I) 
3400 FORMAT(' ',F7.6) 


651 CONTINUE 


GO TO 6000 


4000 WRITE( 30,4010) 

4010 FORMAT(1X,//,T50,'FIXED PHASE RELIABILITY OPTION’ ) 
WRITE( 30,4020) 

4020 FORMAT('-',T38,'PHASE NUMBER' ,T78,' ACTUAL RELIABILITY’ ) 
DO 4030 M=1,NPHASE 
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WRITE( 30,4040) M,AREL(M) 
AOA FORMATC 0 ,141)12,783,58.6) 
4030 CONTINUE 
GO TO 5000 
6000 CONTINUE 


STOP 
END 
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5. SAMPLE OUTPUT FROM DRG FORTRAN Al 


DISCRETE RELIABILITY GROWTH SIMULATION 


MODEL PARAMETER SUMMARY 
NUMBER OF POSSIBLE FAILURE CAUSES 5 


FIXED PHASE RELIABILITY OPTION 


PHASE NUMBER ACTUAL 
1 0 
2 0 
3 0 
4 0 
5 0 
6 0 
7 0 
8 0 
2 0 

10 0 


NUMBER OF PHASES IN THE SIMULATION 10 


PHASE NUMBER NUMBER OF FAILURES IN THE FIRST PHASE 
1 


OVO OnNAMDWNN— 
a eh eh ae eek en oa 


a 


TOTAL NUMBER OF FAILURES 10 


DISCOUNTING PERFORMED USING THE CONSTANT FRACTION METHOD 


FRACTION EACH FAILURE IS DISCOUNTED 0.000 


NUMBER OF TRIALS AFTER A FAILURE BEFORE A DISCOUNT IS APPLIED 


RANDOM NUMBER SEED USED 624712.00 


RELIABILITY 
.398418 
-428109 
-480793 
-539243 
~o0994o 
. 702683 
.798124 
-o999G8 
fy )) Oe) 210 
.990040 


000 


NUMBER OF REPETITIONS OF THE SIMULATION PEREORMED 
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